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ABSTRACT 


This report treats the buckling and_ post- 
buckling behavior of a cylindrical shell that is 
subjected to uniform external pressure p on its 
lateral surface, and an axial compressive force F 
(Fig. A). The force F varies with the pressure p 
m such a way that fF = da2p, in which a is the 
mean radius of the shell and 2 is a dimensionless 
constant. If the shell is immersed in a fluid at con- 
stant pressure p and if the force F results only 
from the pressure p on the ends, \ = z. 

The ends of the shell are assumed to provide 
simple support to the cylindrical wall. Accordingly, 
the radial and circumferential displacement com- 
ponents of the middle surface of the wall vanish 
at the ends. If the ends of the shell are free to 
warp, no other constraint is imposed on the de- 
formation. If the ends of the shell are rigid, the 
axial displacement is constant at either end. Both 
of these cases were investigated. For generality, 
the shell was considered to be reinforced by several 
rings or hoops. 

Only geometrically perfect shells were studied; 
that is, initial dents and out-of-roundness were not 
taken into account. Only shells with a linear stress- 
strain relation were considered. 

If the axial force F is not too great, the shell 
assumes a fluted form when it buckles. This form 
is illustrated by Fig. B, which is a photograph of 
some of Sturm’s test specimens ‘)*. The number of 


ue ed 


Fig. A. Forces Acting on Shell 


Fig. B. Front Views of Buckled Cylinders 


flutes in the buckled form is influenced strongly by 
the ratio L/a, in which L is the length of the shell. 


‘Fig. C illustrates several forms of cross sections of 


cylindrical shells that have been buckled by ex- 
ternal pressure, 

When the axial force F predominates, the 
buckled shell assumes a form in which diamond- 
shaped facets occur (1. Art. 85). This type of 
buckling was not considered in the present study; 
the axial foree / was assumed to be so small that 
the fluted pattern occurs. The admissible range of 
F was not determined, but the fluted pattern usu- 
ally occurs if A does not exceed z. 


* Numbers in parentheses, unless otherwise identified, refer to the 
References at the end of this report. 


Fig. C. End Views of Buckled Cylinders 
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Il. PRELIMINARY CONSIDERATIONS 


1. Introduction 


Experimental data on the collapsing pressures 
of eylindrical shells have been obtained by Fair- 
bairn, Carman, Jasper and Sullivan, Saunders and 
Windenburg, Windenburg and Trilling, Sturm, and 
numerous other investigators. * * ° ® % Theoret- 
ical studies of the problem have been performed by 
Bryan, Southwell, Cook, von Mises, Donnell, 
Sturm, and others. % 1% 122;7) The history of 
these theories (to 1948) is contained in the work 
of Batdorf.”® 

Von Mises and most of the subsequent investi- 
gators implicitly based their analyses on the gen- 
eral principle that a motionless conservative 
mechanical system becomes unstable when the 
value of its total potential energy ceases to be a 
relative minimum. The theory of buckling based on 
this principle is sometimes called the “infinitesimal 
theory” since investigations of relative minima 
require only infinitesimal variations. The buckling 
load determined by the infinitesimal theory has 
been designed by Friedrichs* as the “Euler criti- 
eal load,’ since Euler employed the infinitesimal 
theory in his study of columns. The Euler critical 
loads for elastic cylindrical shells that are subjected 
to external hydrostatic pressure are in close corre- 
lation with experimental data, provided that the 
shells are long in comparison with their diameters. 
However, the Euler critical loads are much too 
high for short thin shells. 

In 1938, von Karman and Tsien‘) called at- 
tention to the fact that an ideal shell can be in a 
state of weak stability, such that a small blow or 
other disturbance causes it to snap into a badly 
deformed shape. Simple examples of this type of 
equilibrium are common. The equilibrium of a coin 
that is balanced on its edge is stable, but such 
weak stability is usually unsuitable for engineering 
design. Similarly, if the center of gravity of a ship 


is so high that the slightest push will cause the 
ship to capsize, the border line of stability has been 
reached. However, this condition has no significance 
for the design of hulls. Analogously, the Euler criti- 
cal load of a shell loses much of its significance 
when snap-through can occur. We merely know 
that the Euler critical load is the upper bound of 
the load that will actually cause failure. 
In this investigation, 

the occurrence of weak p 

stability is manifested by 
the conclusion that the 
pressure required to A 
maintain a buckled form \ 
is frequently much less 
than the Euler critical 
pressure. A pressure- 
deflection curve for an 
ideal elastic cylindrical 
shell that is loaded by 
external pressure has the O 5 
general form shown in 
Fig. 1. The falling part 
of the curve (dotted in 
Fig. 1) represents unstable equilibrium configura- 
tions. Also, the continuation of line OE (dotted) 
represents unstable unbuckled configurations. Actu- 
ally, the shell snaps from some configuration A to 
another configuration B, as indicated by the dashed 
line in Fig, 1. Theoretically, point A coincides with 
the Euler critical pressure E, but initial imperfec- 
tions and accidental disturbances prevent the shell 
from reaching this point. To some extent, point A 
is indeterminate, but it is presumably higher than 
the minimum point C, unless the shell has excessive 
initial dents or lopsidedness. In this report, a 
hypothesis of Tsien“® is used for locating point A. 
The pressure at point C is the minimum pressure 
under which a buckled form can persist. Thus, 


Pressure-Deflection 
Curve 


Fig... 1. 
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if the shell is in a buckled state, and if the ex- 
ternal pressure is gradually relieved, the shell will 
snap back to the unbuckled form when the pressure 
at point C is reached. 

An analysis of the post-buckling behavior of a 
structure determines the buckling load automati- 
cally. For example, an analysis of the form of a 
buckled column reveals that there is no real non- 
zero solution unless the load exceeds a certain 
value, the Euler critical load. Accordingly, in prin- 
ciple, the nonlinear theory of equilibrium obviates 
the need for a special theory of buckling. However, 
as a practical expedient, it 1s usually easier to de- 
termine the Euler buckling load of a structure by 
solving a linear eigenvalue problem than by calcu- 
lating the bifurcation poimt of a curve in configura- 
tion space that represents all equilibrium configura- 
tions. 

Problems of post-buckling behavior of elastic 
shells may be approached mm two different ways. 
On the one hand, we may seek to solve the equi- 
librium equations and the compatibility equations, 
in consistency with given boundary conditions. 
However, in the large-deformation theory of elas- 
ticity, the compatibility equations are an extremely 
complicated set of differential equations, represented 
by the vanishing of a Riemann tensor.@° As Dr. 
C. Lanczos once remarked, “We could not hope to 
solve the general compatibility equations, but for- 
tunately we already know their general solution. 
It is merely an arbitrary displacement vector. We 
should be happy that we know this solution, and 
we should make every possible use of it.” 

When the components of the displacement vec- 
tor are adopted as the dependent variables in a 
shell problem, only the equilibrium equations and 
the boundary conditions remain to be considered. 
The equilibrium equations may be derived by bal- 
ancing forces on a differential element, but, in 
large-deformation theories, the rotations of the ele- 
ments introduce a complexity into this procedure. 
Consequently, the equilibrium equations are ob- 
tained most readily in terms of the initial coordi- 
nates by applying the Euler equations of the 
calculus of variations to the potential energy in- 
tegral. Unfortunately, in most shell problems, the 
equilibrium equations are too complicated to be 
solyed rigorously. Instead of tackling the equilib- 
rium equations directly, we may revert to the po- 
tential energy integral and apply approximation 
methods of the calculus of variations. This pro- 
cedure was employed in this investigation. The 


P- 


theory is accordingly founded on the well-known 
principle that all states of equilibrium — stable and 
unstable — are determined by the stationary values 
of the potential energy. The stable states corre- 
spond to relative minima of potential energy. 

The potential energy 
of the shell is the sum of 
four parts; namely, the 
membrane strain energy, 
the strain energy of bend- 
ing, the strain energy of 
reinforcing rings, and the 
potential energy of ex- 
ternal forces. Articles 3 
to 18, inclusive, are de- 
voted to the derivation of 
the potential energy expression. 

In the development of the theory, the axial, 
circumferential, and radial components of displace- 
ment of the middle surface (wu, v, w) (Fig. 2) are 
approximated by three terms of Fourier series 
(Eq. 11). By using the assumption that the shell 
buckles without incremental hoop strain on the 
middle surface, the Fourier coefficients 14, V2, Va, 
Wo, We, Ws are all expressed as functions of w,. Sub- 
sequently, w, is replaced by a more convenient 
parameter W, defined by W= (n—1/n) w,/a, 
where n is the number of waves in the periphery of 
the buckled shell. It is assumed that W = Wp, cos 
7a/L, where x is an axial coordinate with origin at 
the center section of the shell, and JW, is a constant 
that must eventually be chosen to minimize the 
buckling pressure. The Fourier coefficients wo, Un, 
Us, Us are determined by the calculus of variations 
to minimize the buckling pressure. Accordingly, 
these are finally expressed as functions of W. 


Fig. 2. Rectangular Coordi- 
nates of Shell 


2. Notations 


a = mean radius of the shell 

L = length of the shell 

h = thickness of the shell 

n= afd 

I =moment of inertia of the cross sec- 
tion of a reinforcing ring about its 
centroidal axis 

p = pressure on the lateral surface of the 
shell 

F = axial force that acts on the shell 
(Fig. A) 

\ =a constant, defined by F =da®p 

n = number of complete waves in a cross 
section of the buckled shell 


thr ina 


<Ahew 30; 


EE 
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H = Young’s modulus 
y = Poisson’s ratio 


n ip : 
B= >-V 5 = 0.295804nL/a, if »= 


0.30. 


x = an axial coordinate with origin at the 

center section of the shell 

an angular coordinate (Fig. 2) 

u,v,w = axial, circumferential and radial dis- 
placement components of the middle 
surface due to buckling (Fig. 2) 

V =total potential energy of the shell 
(strain energy plus potential energy 
of external forces) 

AV = increment of potential energy due to 
buckling (Eq. 100) 

U, = strain energy of a reinforcing ring 
(Eq. 99) 

U, = part of the strain energy of the shell 
that results from bending 


> 
ll 


K = constant in the buckling formula, p-> 
=KEh/a 

K; = value of K determined by the infini- 

tesimal theory of buckling 

Ks = value of K determined by the snap- 
through theory of buckling (Tsien’s 
theory) 

., Kis = functions of n and y, defined by 
Eqs. (39), (47), (58), and (67), and 
tabulated in Table 1 

1, 42,01,b2,53,¢1,¢2,¢3 = functions of n, v, and r, de- 

fined by Eqs. (72) and (98), and tabu- 
lated in Tables 2-20 

B,, Bo, Bs = constants defined by Eq. (101) 

W, = A parameter defined by Eq. (36). Wo 
is a measure of the deflection due to 
buckling. 

W = Wo cos 


L 


Kuahka: 


Primes denote derivatives with respect to x. 


Il. POTENTIAL ENERGY OF A SHELL WITH FLEXIBLE ENDS 


3. Membrane Strains 


In this article expressions for the membrane 
strains of the shell in terms of the displacement 
components of the middle surface of the shell are 
derived. 

The shell is referred to rectangular coordinates 
(v,y,2), such that the v-axis is the geometrical 
axis of the cylinder (Fig. 2). The positive x-axis 
in Fig. 2 is directed toward the reader, The circle 
in Fig. 2 represents a cross section of the middle 
surface of the unbuckled shell. The origin of & is 
taken to be the middle section of the shell. 

When the shell buckles, the particle that hes at 
point (x,y,z) on the middle surface is displaced 
to the point (a*,y*,z*). In terms of the axial, 
circumferential, and radial displacements compo- 
nents (uw,v,w) and the angular coordinate @ (Fig. 
2), the coordinates (a*, y*,2*) are given by 


vw=at+u | 
y* =asinéd + vcosé + w sin 6 
z* = acosé ep ey ae 


(1) 


The displacement components (u,v, w) are func- 
tions of x and @. In deriving Eq. (1), we have neg- 
lected the fact that the deformation before 
buckling alters the radius slightly. 

If x and 6 take infinitesimal increments dx and 


dd, the coordinates (x*,y*,z*) take increments 


(dx*, dy*, dz*). These increments are obtained by 
differentiation of Eq. (1); hence 

dx* = (1 + uz) dx + uedd | 

dy* = (vz cos @ + w, sin 6) dx 


+ (acos 6 + v9 cos 6 — vsiné 
+ we sin 6 + w cos @) dé 

dz* = (—v, sin 6 + w; cos @) dx 
+ (—asin 6 — v9 sin é@ — v cos é 
+ we cos @ — wsin @) dé 


(2) 


J 
where subscripts x and 6 denote partial derivates. 

Consider two differential vectors (dx*, dy*, dz*) 
and (8x*, sy*, 8z*), the first being the increments 
of (v*, y*,z2*) when x alone receives an increment 


dx, and the second being the increments of (2*, y*, 
z*) when @ alone receives an increment dé. Setting 
dé = 0 in Eq. (2), we obtain 


0 (ie as 


= (v, cos 6 + w, sin @) dx (2') 
dz* = (—v,sin 6 + w, cos 0) dx 
Setting dx = 0 in Eq. (2), we obtain 
bx* = updé 
dy* = (acosé + v4 cos 6 — v sind 
+ we sin é@ + w cos é@) dé (2”) 


dz* = (—asin# — vg sin é@ — vcosé 


+ we cos @ — wsin @) dé 


The squares of the magnitudes of the vectors 
(dx*, dy*, dz*) and (6x*, dy*, 6z*) are 


(ds*)? = (da*)? + (dy*)? + (de*)? 
(ds*)? = (da*)? + (dy*)? + (62*)? 


Accordingly, Eqs. (2’) and (2”) yield 


(ds*)? = [(1 + uz)? + 022 + w22] (dx)? (3) 
(5s*)? = [we? + a? + ve? + v? + we? + w? 
+2 ave+ 2aw + 2vgw — 2vwe] (de)? (4) 


The initial magnitudes of the vectors (dx*, dy*, dz*) 
and (6x*, 6y*, 6z*) —that is, the magnitudes of the 
line elements before buckling —are 


ds = da, és = adé 
Consequently, 
ds* ri 
( oF ) =14+2u,+uZ2 +02 + w2 (5) 
(2) a1 .2(mte) (2) 
6s 
+ (te y+ (f=. a (6) 


Since the material will not admit large strains, 
the ratios ds*/ds and 6s*/és are approximately 
equal to unity. Therefore, the additive terms in- 
volving, u, v, w on the right sides of Eqs. (5) and 
(6) are small compared to unity. Accordingly, 
ds*/ds and 6s*/és are closely approximated by 


/binomial expansions of the square roots of the 
right sides of Eqs. (5) and (6) in which only terms 
to the second degree are retained. Thus, we obtain 


* 
is ise Us ae ae =e ye 


ds 2 2 
Ose ve + w Ue” 1 ( v — We y a 
is ec a 4 2a? + 2 a Me 


The shell is already strained before it buckles. 
When buckling occurs, line elements in the x and 6 
directions receive incremental strains, Ae, and Ae. 
According to the customary definition of strain, 
these increments are 


: — WS 0s w0S* 208 
Ae, = Abs ) eg — 58 
Consequently, by Eq. (7), 
— 1 2 1 2 
Aéy = Uz + 9 Ve" + 5) Wx 
= Uae Ug? i! ( Vv — We y 
ote a a 2a? nf 2 a (8) 


The shearing strain y.»9 is defined by y.«=cos ¢, 
where ¢ is the angle between the vectors (dx*, dy*, 
dz*) and (6a*, dy*, 62*). Therefore, 


ORO se CUO ap Oe 


‘aaa ds*5s* 
Thus, by Eqs. (2’) and (2”) 
dadé 


ee aaa [(1 + wz) we + (v2 cos 6 + wz sin 6) 


- (a cos 6 + vg cos 6 — vsin 0 + we sin 0 

+ wcos 6) + (—v, sin # + wz cos 6) (—asin 6 

— vp sin @ — v cosé + we cos 0 — wsin 6)| 
Since ds =dzx and 6s =ad6, this equation reduces to 


lie dsianos 
a ere [(1 + we) ue + ave 


+ v, (ve + w) — wz (v — we)| 


Expanding the reciprocals of the right sides of 
Eq. (7) by the binomial theorem, we obtain, to 
first degree terms, 


ds/ds* = 1 — uz, 


6s/ds* = 1 — ee 

a 

Only the first degree terms are needed in these ex- 

pansions, since the first degree terms lead to second 
degree terms in the preceding formula for 7z0. 

Eliminating ds/ds* and 6s/ds*, we obtain, to 


second degree terms 


Us v — We 
Grates = + vz — Wz (2=™ ) 


a 

Ue atv ) 
— Uv, — — | ————_ 9 
2— h(E (9) 
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The axial displacement component w is evi- 
dently small compared to the radial component w. 
Consequently, the term w9?/2a? will be discarded 
from Hq. (8). Also, the terms u,v, and ue(ve+w) /a? 
will be discarded from Eq. (9), since they are 
small compared to the respective additive terms 
vz and us/a. A comparison of the relative magni- 
tudes of v and w is difficult. It has been found 
that the quadratic terms in » exert a predominant 
effect in some problems of buckling of rings. Con- 
sequently, all the quadratic terms in v and w will 
be retained. 

Hqs. (8) and (9) merely give the incremental 
strains due to buckling. The strains just before 
buckling are denoted by e, and e. The initial 
shearing strain is evidently zero. Consequently, 
when the quadratic terms containing wu are neg- 
lected, the complete formulas for the strain com- 
ponents are 


+ 


. 1 1 
ega— €2 + Ur ae iy Ue? ate Dy Wz 


€9 = 9 + une” St >( eel i (10) 
? a 


a 


Wa v — We 
be aa si Bs, = WN (=) 


a 


4. Fourier Analysis of Displacement Components 


Equations (10) express the membrane strains in 
terms of the displacement components of the mid- 
dle surface of the shell. In this article, the displace- 
ment components (1, v, w) of the middle surface are 
expressed in the form of Fourier series in @. Also, 
by the assumption that the shell buckles with zero 
incremental hoop strain, the coefficients in the se- 
ries for the v and w displacement components are 
expressed in terms of a single parameter w. 

In view of the fluted pattern that a buckled 
cylindrical shell adopts, the functions wu, v, w may 
be represented by Fourier series, as follows: 


U = Uo + Wm Cos nO + UW cos 2n0 

+ uz cos 3nO+... 
Vv = 01 Sin nO + w% sin 2nd + v3 sin 8nd +... (11) 
WwW = Wo + wi cos NO + We cos 2n0 

+ w; cos 3n6+... 


Here, n denotes the number of complete waves 
in the periphery. The coefficients w;, v;, w; are func- 
tions of w alone. Only the terms to 3né@ will be 
retained in Eq. (11). 

The membrane strains that accompany buck- 
ling are small, since large membrane strains cause 
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excessive strain energy. This fact is exemplified 
if we deform a piece of sheet metal in our hands. 
Although we can bend it easily, we cannot stretch 
it noticeably. This circumstance implies that the 
middle surface of a buckled cylindrical shell re- 
mains approximately developable, since a wide 
departure from a developable form would require 
large membrane strains. Loosely speaking, the 
“easiest”? way for a shell to buckle is that which 
entails the smallest membrane strains. Conse- 
quently, we introduce the assumption that the 
incremental hoop strain Ae that accompanies buck- 
ling is zero. This assumption does not exactly yield 
minimum strain energy, since the axial strain e, 
and the shearing strain y,» are then too large in 
some regions—particularly the end regions of the 
shell. Consequently, the buckling pressure that is 
obtained with the assumption Ae=O0 is slightly 
too large, both for the infinitesimal theory and 
the snap-through theory. The termination of the 
series in Eq. (11) after the third terms also raises 
the computed buckling pressures, since this ap- 
proximation, like the assumption Aey=0, implies 
artificial constraints on the buckling pattern. 
Eq. (11) yields 


Ve + w Wo 
——— = — + a, cos né + ay cos 2n6 | 
a a 


+ a3 cos 38n0 (12) 
shag = B, sin nO + B. sin 2n6 + Bs sin an 
where 

_ Ins + wi 8s an: 
been a } Lae a eo) 


As was remarked previously, the term w¢? will 
be dropped from Eq. (8). Then Eqs. (8), (12), and 
(13) yield 


Nene ae + a, cos nO + ae cos 2n8 + a3 cos 3nd 
1 : E ; 
+ > (B61 sin nO + Bo sin 2n6 + Bs sin 3n8)? 


With the trigonometric identity, 


sin 7n@ sin jné = z [cos (¢ — j) né — cos (i + 7) n6] 


we obtain, after regrouping terms, 


Wo 


Ae = — 
a 


ote + (BY? + Bo? + Bs*) 
=f > (2a, + BiB2 + BBs) cos nb 
1 


(14) 


+ 7, (20 aoe + Be a oe) cos 2nd 


+ + (2a, — BiB2) cos 3né 


= + (> Bo? + oe) cos 4n 


1 1 
oe BoB3 cos 5n6 — af B3? cos 6n8 


Necessary and sufficient conditions for Ae to — 
vanish are that each coefficient in Eq. (14) vanish. 
Hence, 


B10} Bo = 0, OSs a= 0 
1 1 
a, — {Br = 0, wo + | abr = 
These conditions yield 
i. — eevee Vg = —2nw2, v3 = 0 
n 
— 1 2 2 
Wo & we ee ; (15) 
Wo = 
be 4n? — 1?’ Pee 
or 
en he 5 GS Ware 
a Tip ‘Da (Ane aoe 
V3 = 0 
(n — 1/n)? w? (16) 
Wo } 
4a 
Be (ian a »)2 b) 
Hive! (n — 1/n)? wi ahh 


4a (4n?—1) ’ 

Eq. (16) expresses the coefficients in the v and 
w equations (Hq. 11) in terms of w;. Since the 
curve of a buckled cross section cannot intersect 
itself, the admissible values of w; are restricted to 
a finite range. If w; lies outside of this range, @ 
ceases to be a regular parameter for the buckled 
cross section. 

Eqs. (10), (11), and (16) yield the following 
expressions for the strains (where primes denote 
derivatives with respect to 2): 


é€r = €,% + wy + uw’, cos nO + u's cos 2n6 


+ u’; cos 3né + + (v’; sin n 


+ v's sin 2n8)? + + (wo 
+ w’: cos né + w’s cos 2né6)? 


€ = (17) 


n : 2n x 
1a = sc uy sin né — =e Us Sin 2n8 
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3n ; q 
= one uz sin 3n@ + v’; sin nd 


(17) 


+ v’. sin 2n8 — (w’) + w’s cos nd 


+ Atm ' 


sin 70 


+ w’s cos 2né) ( 
With the trigonometric identities, 


cos n# cos 2nd = (cos né + cos 3né) 


sin n6 sin 2n@ = (cos nO — cos 3né) 


+ 
il 
ON 
: 1 
sin n@ cos 2né = a 


(sin 3n6 — sin né@) 
these equations yield 


1 1 
ce; = [© = wo SF 4 vy’? == 75 Yo? 


/ 
wo? 


thers bert tee] 


= [w: =e = v's + w'ow's 


il 


on 


w' yw! | cos nO 
1 
=i bai +pw? 


_ wie | cos 2né 
(18) 


+ [ws —_ > v’yv’5 + > w! 0's | cos 3n6 


ao [ - + vo? + + we | cos 4né 


n 
Yu = [ - a Ww +01 — w' Bi 
Sry 
+ > W's8, | Sin né 


ae [ - tw + y/5 — Sw! 6: | sin 2n6 


= = w’oB; sin 8nd — on U3 Sin 3nO 


Eqs. (18) are of the form, 
é€x = €, + Co + Ci cos né + C2 cos 2n8 
+ C3 cos 3n6 + Cy cos 4nd 
Yoo = Si sin nO + So sin 2nd + S; sin 36 


| a 


where C and S indicate coefficients of cosine and 
sine terms respectively. 


BEHAVIOR OF CYLINDRICAL SHELLS UNDER EXTERNAL PRESSURE 13 


Eqs. (16) and (18) yield 


1 
Co = wo + a (1 + 1/n?) w'? 


af (2 — 1/n)* werw'? [ An? + 1 ] 
16a? (4n? — 1)? 
(n — 1/n)? [ 3 ] 
G — Hh = — 44 y/,2 
Sart: 4a ar An? Si 
Cp = we + + (1 — 1/n?) w'? 


(n — 1/n)4 ww’? 


ch 4a? (4n? — 1) 
eae, (n — 1/n)? ww’? 20 
et les Aad a Sa 
—(n — 1/n)* wPw'? 
Cn = —— 
16a? (4n? — 1) 
ES ws 
a n 


(n — 1/n)3 
7 4° 


1 Mesa 
[2- ra | wu 


Qn ‘ (n — 1/n) (Qn? +1) wiw's 


So I 2 
% me 2 (An? — 1) a 
5 3n (n — 1/n)? wew'r 
“8s. = oe mn 4 (4n2 — 1) a J 


5. Membrane Energy 


In Sections 3 and 4, expressions for the mem- 
brane strains were derived in the form of Fourier 
series; and the coefficients of the series for the 
displacement components v and w were expressed 
in terms of the single parameter w; We now 
proceed to develop an expression for the increment 
of membrane energy due to buckling in terms of 
the parameter w; and the coefficients wo, wu, Ue, Us 
of the Fourier series for the displacement com- 
ponent wu. 

The membrane energy is (16) 


tan (allot 


+ Qvereg + = ae Pyne | dy (21) 


in which # is Young’s ee vy is Poisson’s 
ratio, and h is the thickness of the shell. Eqs. (19) 
and (21) yield 

L/2 
rine f [2 (€,© =F Co)? St OF as C? ae Ci? 


+ C2 + 2 (69)? + 4y (6, + Co) 9 


Un = 
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= = (he) (SPEC s#) | ie 


in which L is the length of the shell. 
The membrane energy just before buckling is 


aah 
UO = gs ee 


see 2 )2 42 (9)? + 4ye,%e,] da 


This result is obtained by discarding the C’s and 
S’s from the preceding equation. 

The increment of membrane energy due to 
buckling is AUn=Un—Un©. Consequently, by 
the two preceding equations, 


rEHha 
1-7 


200? PO Ce Cee Ce 


AU, = 


L/2 
[ [4 (6. + ves) Co 
0 


= a el — v) (S? +e S.? ae sé) | dx (22) 


The initial axial stress is 


E 
o,0 = ae (e, + veo) 
By statics, 
pores —F eS —pa 
2rah 2rh 


where F is the axial compressive force. We set 
F = dpa? where p is the external pressure on the 
lateral surface. Consequently, 

—)pa 

2rh 


With this relation, the initial strains may be elim- 
inated from Eq. (22). Then Eqs. (20) and (22) 
yield 


ea Mee ee a 


AU, = 


Eh 7 L/2 
a! = i low? + kaw’ ow’? 
1 ae 0 | 

+ ks’ (wi/a)? ws? + kyw's4 


+ ke (w/a)! w'st + ks (wi/a)? w's4 
(Ay) | ko— at? ‘a XP] aw? 


Eh 
tte [> 


l+p 
T 


(23) 


Aap 


+ ks Eh P| (w/a)? w'? 


+ + (1 — v) kg (wi/a)4 w’? 


_ 2dap (1 — 9)? | 


7 ee Be 


+ THe w+ K +0) 


|e 


where ky, ke, .. . 
X, Tf respectively represent the integrals that con- 
tain uw, uw, and us. These integrals are 


L/2 
Wi 
ee if [ws = kgu’; — w'? 
0 a 


eae) ue i) 
2 
~ 5 (=) ko >) wilde (24) 


L/2 WN? 
= [ [ws A oii — a) (=) us? 
0 a 
bs 
+ kiuw'sw'? + kisu’s (=) w' 2 
ws | dx 
L/2 9 n 2 
= 12 a= == Le 2 
if [we+ Sa (2 ws 


W1 
+ kisu’s re 


+ el oa v) kg 


r 
wy? 


2 
— (1 =) ku (SE) wae 28) 
The constants k; are defined by 
pie ues 
; 16n* 
eh nek. OM ee 
e656 (4n? — 1)! [2 (32n4 — 12n? + 3) 
+ 17 (4n? — 1)?} 
Sh tw A/a , , 
kg 16n? (4n? =a 1)? [96n + 44n 
— 17n? + 5] 
1 
ears 
3 ijn ee 
amr Bh loeRe AC UE OZ ry 


lvoe 


"An? (4n? — 1)? 


+ 45n? — 6) 


(n — 1/n)® 
8 (4n? — 1)? 


(60n® — 108n4 


(27) 


kg = (32n4 — 24n? + 5) 


ee 


are functions of n only, and W, 


- nS ot, aaenieedl pe Fp 


_ _ (nm — 1/n)? (8n? + 1) 

3 2 (4n? — 1) (27) 
n(n — 1/n)3 (8n? — 3) 

4 2 (4n? — 1) 

sa + (ay ae 

eae 1 /n)* 

~ 2 Gn? — 1) 
(n? — 1) (2n? + 1) 

ne 4n? — 1 

2 W@W oP 

~ 2(4n? — 1) 

= 3n (n Tr 1/n)? 
Ant — 1) | 


The notations Ky, Ks, ..., are reserved for cer- 
tain combinations of the quantities ky, ke, ete. 


6. Potential Energy of External Forces 


The potential energy of the external forces 
consists of two parts, the potential energy of the 
axial force F and the potential energy of the lateral 
pressure p. If the ends of the shell are rigid, the 
Fourier coefficients uw, wz, us vanish at the ends. 
Then the increment of potential energy of the 
axial force is simply 2F'uo(L/2). 

If the ends of the shell are not rigid, the poten- 
tial energy of the force F depends on the way in 
which the axial load is distributed. We assume 
that it is distributed so that the axial stress o, in 
the cylindrical wall is constant at the ends, «= 
+L/2. Then the increment of potential energy of 
the force F is 


Qa 
AG [ aho,(L/2)u(L/2)d0 
0 


Since o,(L/2) is constant, Fy = —2raho,(L/2). Con- 
sequently 


Qr 
AQe = -f u(L/2)d6 
T Jo 


With Eq. (11), this yields AQy =2Fu(L/2). This 
is the same result that was obtained for a shell 
with rigid ends. 
Since F=)a’p, the preceding formula yields 
AQp = 2a? p09 (L/2) 


Since, by symmetry, wu vanishes at the center sec- 
tion, c=0, this equation may be written in an in- 
tegral form, as follows: 
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L/2 
AQ; = 2na°p [ w' pdx (28) 
0 


To calculate the potential energy of the lateral 
pressure p, we must determine the area A* of a 
cross section of the buckled shell. The intersection 
of the plane «=constant with the middle surface 
of the buckled cylindrical wall is represented para- 
metrically by y* =y*(0), 2* =z*(0). These functions 
are given explicitly by Eq. (1). 

The area enclosed by the curve y*=y*(6), 
2 =2"(@) is 


Qa 0z2* 
At = ait jo ——— 0 29 
Ne aaey. (29) 
The sign on the right side of this equation is nega- 
tive, since the positive sense of @ runs clockwise. 
Eq. (29) is a special consequence of Green’s the- 


orem. 
By Eq. (1), 


y* = asinéd +vcosé + wsiné 
= —asin dé — vsin 6 — v cos 8 
+ we cos 6 — wsin 


With Eqs. (11) and(15), these equations yield 


" : Wi : 
y* = (a + wo) sind — a 08 #sin n@ 
3 2nwo A 
+ w;sin 6 cos né — “Tea cos 6 sin 2né 
Wo : ‘ 
eT sin 6 cos 2n@ 
—dz* : ; 
Cree oi (a + wo) sin 6 — wo sin 6 cos 2né 


+ (n — 1/n) wi cos 6 sin né 
Consequently, if n>2, Eq. (29) yields 


TWo- 


1 a 
1 es 1 (a se Wo)? — — r(1 —1/n*)w? ~ 2(4n?—1) 


2 
By means of Eq. (16), wo may be eliminated from 
this equation. 

The effect of the deformation before buckling 
on the incremental cross-sectional area will be neg- 
lected. Then the increment of cross-sectional area 
due to buckling is AA = A*—7za?. Consequently, 


wy 
AA=qrT [ = fia? + key @ ] (30) 
where 
BESS paw 
kis = 2 (n 1) 


kay = aq = Yn) [1 — sary | a 


Although Eq. (30) has been derived for n>2, it 
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remains valid for n=2, as we see by carrying out 
the integration specifically for n=2. 
The increment of potential] energy of the lateral 
pressure due to buckling is approximately 
L/2 
AQ, = 2p i! AAdx (32) 
Eq. (32) implies the approximation that the axial 
displacement w does not influence the work of the 
lateral pressure p when the shell buckles. 
The total increment of potential energy of the 
external forces due to buckling 1s AQ=AQ;?-+AQ,. 
Consequently, Eqs. (28), (80), and (82) yield 


ale Ge! 


AQ = 2a? rf w'odx + 27a? Pf 


alee 


7. Elimination of uy from the Increment of Total 
Potential Energy 


(33) 


The increment of total potential energy due to 
buckling is AV =AU,,+AQ+AU), in which AU» is 
the increment of membrane energy, AQ is the in- 
cremental potential energy of the external forces, 
and AU; is the incremental strain energy of bend- 
ing. The slight axial bending that exists before 
buckling will be neglected. Accordingly, AU, is 
approximated by the total strain energy of bend- 
ing U,. A detailed analysis of the term U; is devel- 
oped in Section 12. 

Since the term U;, does not depend on the axial 
displacement wu, Eqs. (23) and (33) yield 


wHha rl? 
1 = / Dit oe -+ kau’ yw’? 
0 


= 2 
Dap v) ul] ae 


AV = 


+ ksu’o (w/a)? w'? 


L/2 
+ 2rap [ w'ydx 

0 
+ terms that do not contain w’o. 


By the principle of minimum potential energy, 
the axial displacement w provides a minimum to 
AV. Consequently, w’o minimizes the integrand in 
the preceding equation. A necessary condition for 
the value of the integrand ¢ to be a minimum is 
d¢/du')=0. Furthermore, this condition is suffi- 
cient to insure a relative minimum, since 


ao — 4rKha 
Ow’ 9? 


Salt, 


1-7 


Consequently, 


; ks (wi/a)? ws? — (84) 


1 
4 kaw! 2 


Ht, 


Eliminating w’o from AV by means of Eq. (34), 
we obtain 
aEha 


L/2 
1 = [ IC — Eke) ws! 
elem 
Ages 
2 
ft (4s et kaks ) (+) wt 
4 a 


a xP] ie 
ky Eh | 2? 


—F =») [kt a EP | Coe “ 


1 (Dice Ne <n 
oe r) ks ( ny 


AV =U,+ 


+5 1-9» [he 


nf CS . 
mae ) kee Eh a 
a we WY 4+X+7) 


It will be assumed that 


Tx 
W cos — 
wi L Te 


= w cos 
a n—I1/n L 


(36) 


in which W, is a constant. Observations of buckled 
cylindrical shells suggest that this is a reasonable 
assumption. Since this assumption implies an ar- 
tificial constraint, it possibly raises the computed 
buckling pressure, but it cannot lower it. Eq. (36) 
yields 


L/2 P td 37 tatW 4 
i Ore = "716 (nw — 1/n)A Le 
ie witwiitde = UU 
jot te, ce Dba = Dye ae 
L/2 4/76 [476 
Seek! = bs tas W 0 
I ae 32 (n — 1/n)® L3 
up ew? 
Dea 
0 4(n —1/n)?L (37) 
L/2 274 ]A74 
ek PEL os rat W 0 
' Vi ei earl atte 
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L/2 2 5 6 3 
Ra aD = Ta Wo rv aa 7 
I Se a NF (37) 81 +7) (nm? — 1? 2) 
Ne vLW mk 
da — eae 7 
ee aod (n — 1/n)? 32 (1 +p) (n — 1/n)4 
O odin meee Wir de & = ua [ = 
16 (nw — 1/n)* J 32 (1+ 7) L (4n? — 1)? 
Consequently, by Eq. (35), 4 2n? — 1 ie 3n? — 4 ] 
: - 2 (n? — 1) (4n? — 1) — 4 (nm? — 1)? 
AV =v,+[K.—2* — x, 2 
b 6 ie 4 Ts K wks 4 7 
pe = 64 (1 +7) (n—1/n)® 5121 +2) 
— Kwa®pl | We + Es L 1 9 
= D) 
Hath ae [ (4n? — 1)? An? — I = 2| 
ee Ga a 74 (38 
—K, =" —K, eee pL | Wot (38) i eee) flea 
2a ae 2G 
E ae) kK. Eah ] We (n /n) (n ) 
=e =e 0 Ami) K Az 3rky; 
ae: ‘ atl " 8 (n — 1/n)* 
a hc 
in which 128 2 (An? = 1) J 
(x 1 ks :) These constants have been tabulated (Table 1). 
m \ Ky — ~5- Kg 
K, = — 8 8. Elimination of u, from the Increment of Total 
16 (1 — v?) (n — 1/n)4 Potential Energy 
a 3n° Since the origin of x is the mid-section of the 
258 (Dia) (ne — 1)? Shell, the displacement component wu is an odd 
1 function of x. Therefore, w is an odd function of 
3m? @ rear hs) x. By the principle of minimum potential energy, 
Ky = 256 (1 — »*) (n — 1/ny8 this function provides a minimum to W. When w: 
is eliminated by means of Eq. (86) the equation 
ne 51° for WY becomes 
65536 (1 — v*) (4n? — 1)? ot ao 
ie i: [ws ese Sa uw’; cos —— TE sin? (22 *) 
5 (i — + hiks) : - 
ie a : ne EE 5 oe es ez (eyecare 
= 32 (1 — »*) (n — 1/n)8 Ss tig Sa 2 B) 7 ea 
ae x 1 
5n* in? — 1 sin ~~ cos? ( m) +—( u | dx 
(4n? — 1)? + An? — 1 as 2n | # e - 
aa FA 
256\(1 =p) (ne = 1) where, for brevity, ee 
ms mks De (nell) te ' : 
ea Pare ee OO: With the trigonometric ace : 
se Eels aE sin? (22 ) = : (cos fe sce re) 
oe 32 S = hin)? Be a Tx 37x 
sin TE eos? ( 2 q + (sin — + sin L ), 
3 4n? + 1 | 
= oe ea (4n2 — 1)? (39) this yields 
L/2 2 2,3 
by ke = re gent oo ws( Rube Sat ) 
K; = Sa + = =n Vv [ [« 1 ko Fa wis \ COs > COS 
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peas, = ar ~ (sin 2 


oat ) eines nee). 


ce integration by parts now yields 


een 7 UW sin —— 


+ sin 
Since w(0) 


L/2 2A rx 
— 12 25-2 11) 
v= f [ wit + keut + iL uy sin =| 
a = uy sin ort | ae (40) 
in which 
1 3 2 
Ang [Fame Pete 
—+ (1-9) x 
Gael) 
1 31 2 
p~3[F0- nit ee 
l—vyp 
is 2 J 


For integrals of this particular type, Euler’s 
equation of the calculus of variations is both 
necessary and sufficient for a relative minimum.“? 
EKuler’s equation for WV is 
WNT Sas Gites B 3re 


7 1 4 D sin L (42) 


The general odd solution of Eq. (42) is 


uw", — Puy = 


: AUT eG ee eer Tra 
u.=CL sinh kx @ S| B sin (43) 


in which C is an arbitrary constant of integration, 
and 

a=r+hL?, B= 97? + kL? (44) 
a ieee Eq. (43) into Eq. (40), we obtain 


ACL BL 
=) 27 2 ee Fa ee 
WV = ke L? sinh kL a 4B 
If there are no constraints at the ends that 
affect the axial displacement u, the constant C 
must be chosen to minimize WV. Obviously, this 
condition implies that C =0. Then, 
A[]L BL 
BE dae th) Age, Se 
Evidently, VY is negative and it reduces AV. This 
circumstance might have been anticipated, since 
the introduction of 2 effectively gives the shell 
additional degrees of freedom. 
Let us set 
Kx tk ie mn (8n? + 1) 
ni 8(n — 1/n)3 16 (n? — 1) (4n? — 


(45) 


D | (47) 


meet TF (1 = y) ko 
Kae 16 (n — 1/n)8 
a (1 — v) n (8n? — 8) 
32 (4n? — 1) 
2 Ws vy) 7m eo 
Ku = 2G =e r= a/L 


Then, by Eq. (41), 


A= (Kio —— Kor?) We = KuWo 
= (Kio + 3Kor*) Wo 


il ye iE 
pee - oe 


Accordingly, Eq. (46) yields 


T 


ae = Wen) + Le© + rf(@)] Wet 
+ [falé) — fale) + r'fol€)] Wo (49) 


where 

—7khy? —1 Ky? il 
OS a = es fi(é) = 7 ¢ ee a (++ +5) 
AO=sq aay HO= ios a 


+e —TrKy Ky ee 5 (+ >) 
9. Elimination of uz from the Increment of Total 


Potential Energy 


Like wi, uw. is an odd function of x. This fune- 
tion must be chosen to minimize X. When wy; is 
eliminated by means of Eq. (36), the equation for 
X becomes 


24 fe ty br Gin yy n() 
ae a 7 ([40- L/n®) w a zat ) 
Se ee se t G- Ws Pe 


An? — 
a eee 1) one +), 
2L (4n? — 1) 

- 5 ome 

sin } dx (51) 
Integration by parts yields 

2A ee 
= [ws + Ake us = iF Us. SIN = 


— 28 wz sin L Jax tye (1 
— 1/n?) wu, (L/2) (52) 
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where 
Repel ES — v) (nr? — 1) (2n? + 1) 
4 4n? — 1 
+ xr? (1 = 1/nt) | o? 
(53) 
ba ue (n — 1/n)4 
eS 4n? — 1 a 
= P= OLE. 


Kuler’s equation for the integral X is 
A Qrx B Arx 
‘gga 2 = : —. s 
ws — Aku. 7, sin] 7, sin | 


_ The general odd solution of this equation is 


‘ AL Qrz , BL . 4rax 
U2 =CL sinh 2kx+ 4 iB + Te sin 5 (54) 
where 
= 7+ kL, y = 40? + k2L? (55) 
Substituting Eq. (54) into Eq. (52), we obtain 
A’L BL 
& 2012 ot a ns 
X = kL?C? sinh 2kL iGo 16y 
— a sinh kL + Lands sinh kL 
oe — 1/n?) w&C sinh kL (56) 


3L 
If the ends are not constrained against warping, 
C must be chosen to minimize X. Then Eq. (56) 
yields 


De ae ee Ba Tee 
Sages (6a 164 ase Gin) a 
2 
ee ae are] tanh kL (57) 
a ¥ 
Set 
pe (1 — v) n? (2n? + 1) 
a 4 (nr? — 1) (4n? — 1) (58) 
T° pS T° , 
ea Cae ey? MS avant | 


Then, by Eq. (53), 
ik, 2 Kee Be Koes! 


(59) 
1 n Lae 
Cg en aed | 


Consequently, to sixth degree terms in Wo, 


oe Witeilver Wildl.’ 860) 


1-v L 
where 


gilt) = 


Tv (Kip a Ky3r?)? 
160 —v)a 


| (61) 
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7 tanh 2& | ar 
161 — »)é— L 2 (mn? — 1) 


a Ki ale Ky3r? ] 


a 


_ 


(61) 


m Kir? tanh 2 [ ar 
4 (1 — v?) yé 2 (n? — 1) 


aX Ky + Kas’? ] 


a 


go(€) = 


Numerical calculations show that the eighth de- 
gree terms in W, are entirely negligible. 


10. Elimination of u; from the Increment of Total 
Potential Energy 


The Fourier coefficient w3 is an odd function 
of x that must be chosen to minimize T. When w; 
is eliminated by means of Kq. (36), the equation 
for T becomes 


ak 402 ow? TX TL 


st — ie | w 3° + 9k us eg 5 ee u's sin? — cos — 
L L 


+(1- y) hae ~ us cos? = sin 2 ] da 


With the trigonometric identities, 


sin? ( TY eos LES (cos z= C08 7) 
L ee © 5 Sy ai 

5 ( rx ) sin 7% ¢ (sin NOD ig SUE ) 
cos’ | > — ) sin = | et 


this equation yields 


L/2 2) De 
i i [ ws! + 9k?u3 + aes if 


0 


OTe ahy50* ( Sete 
= HO ant) + (1 — ») az, v3\sin=— 
+ sin ane y] dx 


Integration by parts now yields 


L/2 2A se hat 
TP = i [ ws! + 9k?u3? + L U3 Se 


2B | Ome 
7 vssin | ] dx (62) 
where 


3 
A= a [wksar? + (1 — v) kasl 


(63) 


3 
B= [Sa2kur? — (1 — v) kis], r = a/L 


8 


Euler’s equation for the integral T is 


uz’ —9 k?uz = - sin L z sin 7 
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The general odd solution of this equation is 


; AL . 1x2 
us = CL sinh 3kx — ae sin a 
BG Oo TL 
+9, sin | (64) 


where C is a constant of integration, and 


a=r+RP, 6=r+9rL (65) 
Substituting Eq. (64) into Eq. (62), we obtain 
Ser ; A7L  S BL 
T= ig kL? sinh 3kL 46 36a 
To minimize T, we set C=0. Then, 
APL BL 
ERR Teg oe (68) 
The following notations are introduced: 
K x rky4 ren | 
1 Ee Weiee Pan Aaa? 
8 (n — 1/n) 16 (n 1) (An 1) (67) 
peas mil J) his om (le yin 
Pe Sn 1 n)* 5, 32 Ae a) 
Then, by Eq. (66), 
T ene aint 
where 
e Tv (Kya se Ky;)? 
g(&) = al 5 
at gag OP ml NO 
+ (3Kyar Kis) ] (69) 
9a 
— (on = ftp 
Obes 7: 2 


11. Further Simplification of AV 


In Sections 5 to 10, we have developed ex- 
pressions for the various components of the in- 
crement of the total potential energy due to buck- 
ling. In this article, we proceed to further simplify 
the expression for AV (Eq. 88). 

Set 


EG KG he Bee. Ue" n Le? oa 
p= KE a fs ay by = AJ 5 (70) 


Then Eqs. (88), (49), (60), and (68) yield 


AV) a ° 
Fahl a 


+ (bo + Kaz) Wot + b3Wo® (71) 


in which 
a, = Kar? + Ke 
2 = — K;\r? + Ki; (72) 
= fi + Ker? 


a 
S 
| 


se 
| 


be = fo — (Kz — fs) r? + Kart — gi (72) 
bs = fa + (Ka — fs) 1? 
+ (Ks + fe) r*# — 2 —@ 


An eighth degree term in Wp» has not been included 
in Eq. (71), since numerical calculations show 
that it is very small. The constants (a1, a2) come 
from the potential energy of the external forces, 
and the constants (b1, be, bs) come from the mem- 
brane energy. 

The constants Ki, Ke, ... have been tabulated 
for y = 0.30. (Table 1). With these data, the con- 
stants a1, 2, bi, be, bs have been tabulated for \ = 7. 
(Tables 2 to 20). It is to be noted that affects 
only a and ay; it does not enter into bi, be, bs. In 
the construction of the tables, it was convenient 
to treat € as an independent variable. The corre- 
sponding values of L/a have been included in the 
tables. 


12. Strain Energy of Bending of an Elastic Cylinder 


In Eq. (71), AV is expressed in terms of the 
single parameter W, and the strain energy U; due 
to bending. By the theory 
of this article, U; is also 
expressed in terms of Wo. 
Also, the strain energy of 
reinforcing rings is ex- 
pressed in terms of Wp. 
Thus, the increment of the 
total potential energy is 
reduced to a function of 
the single parameter Jp. 

For an elastic cylindrical shell, the strain en- 
ergy of bending, per unit area of the middle sur- 
face is @® 


Fig. 3. Arc of Buckled 
Cylindrical Shell 


+ D [x22 + Ke? + 2vK.K9 + 2 (1 — v) 7] (78) | 


where x, is the change of curvature in the lon- 
gitudinal direction, kp is the change of curvature in 
the circumferential direction, and + is the local 
twist. The flexural rigidity D is defined by 


Eh’ 

es, 12 (1 — »?) tee 

To express x» in terms of the displacement 
components (v,w) of the middle surface, we refer 
to Fig. 3. The are in the figure represents a part 
of the cross section of the middle surface of a buck- 
led cylindrical shell. Since, by assumption, there 
is no incremental hoop strain due to buckling, 
ds* = ds = adé. Also, by Fig. 3, Rdg@ = ds, where 


FR is the radius of curvature of the are. Further- 

more, we see by Fig. 3, 
sin ¢@ = — de*/ds, cos¢ = dy*/ds_ — (75) 

Differentiating these equations, we obtain 


at do dPy* 


a , a= 


ds ds? ? ihre URE i) 
By geometry, 1/R = d¢/ds. Consequently, 
i SG Cie Gis (77) 
R cos ¢ sin ¢ 
Introducing the parameter @ we obtain by Eq. (77) 
Le Le Se ie dee 
R a Yor ih eye 
Since (y*)? + (z*)? = a’, this equation yields 
1 26*Yon* — Yo*Ze0* i 
Ree as ie) 
By Eq. (1), 
Ye" = acosé + v9 cos# — vsin# 
+ we sin é + w cos é 
zo* = — asind — vgsinéd — vcosé 
+ ws cos? — wsin 8d (79) 
Yoo® = —asind — Qvasin 6 + vo9 cos 6—v cos 0 
+ 2 we cos 6 + wee sin @ — wsiné 
Zeo* = —acosé — 2v9 cos @ — veg Sin O-+0 sin 8 
— 2we sin 6 + wes cos é — wcosé 


Eqs. (78) and (79) yield 


Lge = \[a + (vp + w)] [a + 2 (vs + w) 


R 

— (woe + w)] — (v — we) [(ree + we) 

— (» —w)]} (80) 
By Eqs. (12) and (16), 


Ve + W = Wo — Wo COS 2N0 
v — We = (n — 1/n) wi sin nO 
Wee + wW = wo — (n? — 1) wi cos né (81) 
— Wo cos 2nd 
Veg + We = Znwy sin 2n6 
1 


oa Baa) (e0)-and (61) yield 


Since x= 
Its a’ 


hey = |[ Zar» — 3 Wor + + (n — 1/n)? wi | 
+n(n — 1/n) (a — 5 ro) wi cos né 
= [2 (a + wo) wo + + (n — 1/n)? wi | 


cos 2né + + n(n — 1/n) wows cos 3n8 
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ae > Wo” COS 4n6} 
With Eq. (16), this yields 
ake = nW cos nd + + nW?* sin?né cos n6é 
+ ma W* sin? n6 (82) 
where W = (n — 1/n) wi/a 
Eq. (82) yields 


ee en WE See [ 72 1 tas 
9 Da [rerao = FL =e nw eg eee 
1 271776 35 vs | 
+ 39°" ws + i024 W (83) 


Subsequent calculations show that rarely, if ever, 


does W exceed the value ah Consequently, the 
eighth degree term in W is quite negligible, and it 
will not be included in the subsequent analysis. 
Eq. (86) may be written 
W = Wo cos + (84) 
Eqs. (83) and (84) yield 


1 cane ee oh heh [ ae 
5 Da [i . i uid) = eH | 


DOL ie OS ‘| 
= 16 Wott 256 Wo° | (85) 


Eq. (85) represents the principal part of the 
strain energy of bending. Since the longitudinal 
curvature is a minor effect, we shall use the ap- 
proximation, 


Ky = — Wer (86) 
Then, by Eq. (11), 
Kr = —(wo” + wi” cos né + we” cos 2nd) (87) 
Consequently, 
rHha ie 


ih Qa 
= Daf kx = (Qw"o? + wi? + ws) 
0 


2 24 (1 — »*) 
With Hq. (16), this yields 
i | Ww” 
24 (1 — v) |(n — 1/n)? 
1 1 
"5 4. [2 “s (4n? — ae 
[W'+2WW?w"+ ww) (88) 
Eqs. (84) and (88) yield 
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1 FE 
= (4n? — 1)? iW of ee) 
Eqs. (84) and (87) yield 
le "\Wi+w’ (nw 


Qrv 
a 


[ 2excnado = — 
0 


i yi zw vs | 
+e nW g weW 


With Eq. (16), this yields 


Qe 
[20eeKsd0 = —2e» l= Ww? + Ww") ws 
0 


Ww” 1 ~ 
oF n—1/n (nw si ae 
We P 7/2 vu 
Seer )| 0) 


Eqs. (84) and (90) yield 
wvEhia J n> 
24(1—»?)L | n?-1 
3n? ie 
a 32(n?—1) a Ga [3 


err] 4} oD 


The twist 7 is defined by +r = 0¢/dx (Fig. 3). 


if L/2 
~ Da Te. dx ie Qe kd = We 


Since sin ¢ = — dz*/ds and cos ¢ = dy*/ds, (see 
Kq. 75), 
TCOS¢ = ae ee 
dxds Ged 
: 07y* 1 
rTsin¢@ = — oe aie Yo 
Hence, 
T COs? ¢ = —s Yorer9*, TSin?d = a 2e*Y 26" 
a x8 > ae x 
Adding these equations, we get 
ar = Ze* Ysa" — Yo*Zz0* (92) 


Eqs. (1) and (92) yield 


a’r = a(v’ — we’) + (ve + w) (v’ — wy’) 
— (v9 + w’) (v — we) (93) 


Consequently, by Eqs. (16) and (82), 


ar = (« W' 3. Ww’) + 2 Wwe) sin n6 


. > (Wale aaW ane (94) 
Hence, 
25 i 
2a (1 — ») [70 = mot 
0 a 
3 
=e: Ww" 0 + 3 Ww), 
1 yf al 
+ (Wu's = ww"*| (95) 


By Eq. (16), wo = -+ aW?. Consequently, 
Qa 

9a (Cl ») f 72d9 = 2x (1 —») a(we 
0 


+2 wwe + 
Hence, by Eq. (84), 


ee rr) 
aa: aM 


LJ2 a pee 
2a(1—») [| def ap eee mG 
+5 Wilt525 W'] (96) 


By Eqs. (73), (85), (89), (91), and (96), the 
total strain energy of bending of an elastic cylin- 
drical shell, excluding the strain energy of rein- 
forcing rings, is determined by 


a =F a Wo? (er + c2Wo? + c3Wo') (97) 
where path 
2 See es es Qp2 
ae wd = P) [nt + 2am8(1 + 5 ee 
rin? 4 
+a 


T 


a= Sqm [wt (1+ az) (98) 


ae our (2+ | 
Tv 


2 27.2 fae 
12288 (1 — [5 n? + 27°77 (5 17p 


Sy 
a An? — 1 )] 
The constants ¢1, ¢, c; have been tabulated (Ta- 


bles 2 to 20). 
The strain energy of a reinforcing ring is 


en 


1 Qe F 
i= zo ET xg2d6 


in which J is the moment of inertia of the cross 
section of the ring. Consequently, if J is constant 
for a ring, Eq. (83) yields 


se 


Ue mW? + eW! + ane] (99) 
The eighth degree term in W has been discarded 
from Eq. (83). 

Eq. (84) determines the value of W for any 
ring. The strain energy of all reinforcing rings is 
represented by }\U,, where the sum extends over 


all rings. 


It is assumed in the derivation of Eq. (99) that 
the centroidal axis of the ring coincides with the 
middle surface of the shell. Actually, this condition 
rarely occurs in practice. Possibly Eq. (99) can be 
retained for off-center rings if EJ is modified to 
‘take account of an effective width of shell that 
acts with the ring. However, this problem is not 
examined in the present investigation. 
Eqs. (71) and (97) yield 


me AV 
| Eahl 


= (B, — Ka) Wo? + (Bs + Kas) Wot 
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ay eee (100) 
Bee hl 
in which 

he h? 

By = by + t1—,-, Bo = by oP Case, 
a a 
h? 

B, = bg + C3 a (101) 


The constants a1, dz, bi, be, bs, C1, 2, ¢3 have been 
tabulated (Tables 2 to 20). The last term in Eq. 
(100) representing the effect of reinforcing rings, 
is a sixth degree polynomial in Wo, in which only 
even powers occur (see Eq. 99). 


Ill. POTENTIAL ENERGY OF A SHELL WITH RIGID ENDS 


13. Shell with Rigid Ends 


In the preceding analysis, it has been assumed 
that the ends of the shell are free to warp out of 
their planes; that is, that the ends of the shell im- 
pose no restrictions on the axial displacement w. 
If the ends of thé shell are rigid plates, the Four- 
ier coefficients i, wm, us (See Eq. 11) vanish at 
the ends. Accordingly, the functions VY, X, and T 
must be modified. Numerical computations show 
that the effect of ws is quite small, and conse- 
quently the function YT will be discarded. 

Eqs. (40), (41), (42), (48), (44), and (45) re- 
main valid. However, the constant C must be 


chosen so that uw; vanishes for x = L/2. Conse- 
quently, 

8 gel B kL 

C B esch 5 (102) 
Eqs. (45) and (102) yield 

LA? LB? 
Es 4a 48 
45 (4. = =) Le cothé (103) 
a B 


in which € is defined by Eq. (48). Consequently, 
W 


>, = A’F: (6) — 2A BPi(E) + BPPs(é) (104) 
where 
1 
2 coth & — (2 fe + 7) 
ples Ge? + 4EP 
2— cothé (105) 
F(t) = (a a AE) One aes) 
; _ 2écothé — (+ 9 1?/4) 
es (Om + 42) 
Hence, 
ES = W'yal 
ey Veil [y»(é) 
+ 1°ys()] Wot + [Ya(é) 
= PVE) po(£)| Wo (106) 
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where x 
KiFy 
W(é) = a 
ae. 27 KyKu (FP, aa F) 
Yo(E) ihe yp 
_ nKoKu (Fi + 3F) 
W3(E) 1—,r 
(107) 
ie aK? (Fy — 2F. + Fs) 
W(€) 1 — yp 
2 Qn Ky y= 2s — 8p 
Ws(E) Te D 
= aK? (Fy + 6F2 + 9F3) 
We(E) ie pe 


It may be shown that the expression on the 
right side of Eq. (104) is a negative definite quad- 
ratic form in A and B. Consequently, VW always 
reduces the membrane energy. For a shell with 
rigid ends, the functions Wi, ... Ys replace the 
functions fi, ... fs» of Eqs. (50) and (72). 

Turning attention to the function X, we ob- 
serve that Eqs. (54) and (56) remain valid. The 


end condition, w = 0, obviously requires that 
C = 0. Consequently, Eq. (56) yields 
2 2 
g= Se (108) 
Hence, is xf 
x 
fa (109) 
where 
ee ee (110) 


16 (1 — v*) (re? + 42)? 
Accordingly, the function X, replaces the function 
¢, of Eqs. (60) and (72). The functions ¢. and g 
are discarded from Eq. (72). The second term in 
Eq. (108) has been neglected, since B* is of eighth 
degree Wo. 

With the above modifications, the preceding 
theory applies for a shell whose ends are rigid 
plates. 


IV. PRESSURE-DEFLECTION RELATIONS 


14. Load-Deflection Curves 


For a given shell and a given value of the pres- 
sure p, we may plot a graph of AV/(HahL) versus 
W, by means of Eq. (100). The forms of the 
graphs, corresponding to several values of p, are 
illustrated by Fig. 4. The pressures indicated on 
the curves are such that p, < po <ps3 < ps. The 
minima on the curves represent configurations of 
stable equilibrium, and the maxima represent con- 
figurations of unstable equilibrium. If p < p,, the 
unbuckled state is stable, since the configuration 
W, =0 then provides a relative minimum to the 
potential energy. However, if p 2p,, the unbuckled 
state becomes a configuration of maximum poten- 
tial energy; hence, it is unstable. Accordingly, p,4 
is the Euler critical pressure. 

We may pick the maximum and minimum 
points from the curves of Fig. 4, and thus plot p 
versus W,. The resulting curve, illustrated by Fig. 
5, represents all equilibrium configurations. Fig. 5 
effectively a load-deflection curve for the 
buckled cylinder, since W, is roughly proportional 
to the incremental deflection at the center of a 
lobe due to buckling. The intercept of the curve 
with the p-axis is the Euler critical pressure. The 
falling part of the curve (dotted in Fig. 5) rep- 
resents unstable configurations, since the points 
on this part of the curve correspond to maxima 
on the curves of Fig. 4. The rising part of the 
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Fig. 4. Increment of Potential Energy 
versus Deflection Parameter 
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curve represents stable configurations, as the points 
on this part of the curve correspond to minima on 
the curves of Fig. 4. The minimum ordinate on the 
curve of Fig. 5 is the lowest pressure for which 


p 
= 
‘N bee 
regu! TP See 
~ 
Euler SS 
SS 


critical pressure 
W, 


Fig. 5 (above), Pressure-Deflec- 
tion Curve for Ideal Shell 


Fig. 6 (right). Intersecting 


2 W, 
Pressure-Deflection Curves 2 


the shell will not snap back if it is forced into the 
buckled form. It is equal to p,, if the curve corre- 
sponding to p, in Fig. 4 is considered to have an 
inflection point with a horizontal tangent. 

We may plot Fig. 5 directly by means of Eq. 
(100). The poits on Fig. 5 are solutions of the 
equation dy/dt=0, where, for brevity, 
y = AV/(EahL) and « = W.?. If there are no re- 
inforeing rings, this condition yields 


— Bi + 2B.Wo? + 3B3Wo! 
% Ot = 2a. W ¢? 

In view of Eq. (99), the effect of reimforcing rings 
is merely to modify the coefficients 6,, bs, bs. Con- 
sequently, the form of Eq. (111) remains valid for 
a cylinder that is reinforced by elastic rings. Fig. 5 
is a graph of Eq. (111). It is irrelevant whether 
the ordinate is p or K. 

A curve of the type illustrated in Fig. 5 corre- 
sponds to each value of the integer n. It is neces- 
sary to choose n by trial to provide the minimum 
buckling pressure. In some cases, the curves corre- 
sponding to two consecutive values of n intersect, 
as illustrated by Fig. 6. Then the number of lobes 
in the final buckled form may possibly be different 
from the number of lobes in the infinitesimal pat- 
tern that precipitates buckling. However, since the 


K 


VD = Koby ce. Cai) 
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collapse of an ideal shell is a sudden phenomenon 
that carries the shell over the unstable region onto 
the rising part of an equilibrium curve (Fig. 6), 
dynamical processes undoubtedly play an impor- 
tant part in determining the final pattern. 


15. Tsien Critical Pressure 


The curve corresponding to p; (Fig. 4) 1s con- 
sidered to have an inflection point with a horizontal 
tangent. Then the curve corresponding to any value 
of p in the range p, < p < ps possesses two rela- 
tive minima, one representing the unbuckled state, 
and the other representing a buckled form. Conse- 
quently, if p > p,, the shell will maintain a buckled 
form if it is initially forced into that condition by 
external disturbances. Since initial disturbances 
and imperfections always exist, von Karman and 
Tsien originally conjectured that p, is the maxi- 
mum safe pressure. 

Tsien“ later concluded that, although p, is the 
greatest lower bound for the pressures at which 
buckled configurations can persist, there is little 
danger of a shell passing into a buckled configura- 
tion unless, in doing so, it loses potential energy. 
In Tsien’s words, “The most probable equilibrium 
state is the state with the lowest potential energy. 
— This principle of lowest energy level is verified 
by comparing experimental data with theoretical 
predictions. However, in view of the prerequisite 
that arbitrary disturbances of finite magnitude 
have to exist, the buck- 
ling load determined by 
this principle may be 
called the ‘lower buckling 
load.’ The classic buck- 
ling load that assumes 
only the existence of dis- 
turbances of infinitesimal 
magnitude may be called 
the ‘upper buckling load.’ 
Of course, by extreme 
care In avoiding all dis- 
turbances during a test, 
the upper buckling load 
can be approached. The 
lower buckling load, however, has to be used as a 
correct basis for design.” 

According to Tsien’s reasoning, the pressure De 
(Fig. 4) is the maximum safe pressure. This is the 
pressure at which the potential energy of the un- 
buckled form ceases to be an absolute minimum. 
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Fig. 7. Pressure-Deflection 
Curve for Imperfect Shell 
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Fig. 8. Statistical Distribution 
Curve for Imperfect Shell 


Fig. 9. Determination of Tsien 
Critical Pressure 


It will be designated as the “Tsien critical pres- 
sure.” The curve in Fig. 4 that corresponds to pz 
is tangent to the axis of W, at a point to the right 
of the origin. 

The buckling pressure of an imperfect shell 
poses a statistical problem. Load-deflection curves 
for imperfect shells have the general form shown 
by Fig. 7. This figure is to be contrasted with the 
load-deflection curve for an ideal shell (Fig. 5). 
Donnell has emphasized that the designer is con- 
cerned principally with the maximum pressure on 
the load-deflection curve (denoted by p; on Fig. 
7). Since the falling part of the curve (dotted in 
Fig. 7) represents unstable equilibrium configura- 
tions, the maximum point lies at the boundary of 
the stable range. Therefore, p; is the Euler critical 
pressure for the imperfect shell. This pressure may 
be expressed conveniently as a fraction f of the 
Euler critical pressure ps for a perfect shell; that 
is, f = p;/ps. Since p; depends on initial imperfec- 
tions in the shell, tests of a large number of shells 
with the same dimensions would lead to a statistical 
distribution curve of the general form shown in 
Fig. 8. The ordinate ¢ of this curve is defined by 
the condition that ¢df is the probability that a 
random shell will fall in the interval (f, f + df). 

The specification of a safe pressure is somewhat 
arbitrary. Under some circumstances, an operating 
pressure would be considered safe if 95% of all 
specimens would fail above that pressure. In other 
cases, the safety limit might be raised to 99%, or 
some other value. Tsien implied that his definition 
of the lower critical pressure provides a value that 
lies near the maximum safe pressure (Fig. 8). At 
present, this conclusion is largely conjectural, but 
since Tsien’s eritical pressure affords a ready em- 
pirical criterion for safe design of shells, it has 
been charted by Euler critical pressure (Fig. 13). 

The Tsien critical pressure is determined by the 
equation AV = 0. If we select the intercepts of the 
curves of Fig. 4 with the wpo-axis, we can plot 
the resulting relation between p and Wp». The graph 
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has the general form shown in Fig. 9. The mini- 
mum ordinate of the curve is the Tsien critical pres- 
sure, and the corresponding value of W, determines 
- the deformation of the buckled shell, if the applied 
pressure equals the Tsien critical pressure. The 
intercept of the curve with the p-axis is the Euler 
critical pressure. Although Fig. 9 looks like Fig. 5, 
the two curves are distinct, since they are derived 
by different formulas. Fig. 9 is not a graph of 
equilibrium configurations; it merely serves to show 
how the Tsien critical pressure may be computed. 

If there are no reinforcing rings, and if the shell 
is elastic, the equation AV = 0 yields 


By + BW? + BsWo' 


1 — We 


K= 


(112) 


The form of Eq. (112) remains valid for a cylinder 
that is reinforced by elastic rings, if the coefficients 
bi, bz, bs are modified suitably (See Eq. 99). 

Plotting K versus W,) by means of Eq. (112), 
we obtain a curve that is essentially equivalent to 
Fig. 9, although the ordinate is AK rather than p. 
The intercept of the curve with the K-axis is de- 
termined by setting W,=0 in Hq. (112). Conse- 
quently, the value of K corresponding to the Euler 
critical pressure is 

K; = By/ay (113) 

The notation AK; denotes the value of AK that is 
obtained by the infinitesimal theory of buckling. 

The value of K that corresponds to the Tsien 
critical pressure, denoted by K,, (the subscripts 
“st” denote snap-through) is the minimum value 
of the function defined by Eq. (112). To minimize 
K, W,? must be a root of the equation, 


1 BE a2By ] a2,W o* 
rs 2 B; a,B3 2 2a4 
Although Eq. (114) may be solved by the quad- 
ratic formula, it is usually poorly conditioned for 
this type of solution, and it is most easily solved 
by iteration. The procedure is to obtain an ap- 
proximation of W,? by neglecting the fourth degree 
term on the right side of Eq. (112), and to use 
this approximation to refine the first approxima- 
tion. The process of refinement may be iterated, 
and it converges quite rapidly. The value of W,’, 
determined by Eq. (114), must be substituted into 
Eq. (112). Thus, the Tsien critical pressure is de- 
termined. It is represented by p = K,,Lh/a, where 
the subscripts “st”? denote “snap-through.” The 
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Tsien critical pressure aN 


Fig. 10. Graphs of f(p) and $(p) 


buckling coefficient A: is plotted versus W, in Fig. 
13 for values of a/h of 100 and 1000 and A=7z. 


16. Effect of Assumptions on the Tsien Critical 

Pressure 

It is well known that arbitrary assumptions 
about the deformation of a structure cause the 
computed value of the Euler critical load to be 
too high. The same conclusion applies for the Tsien 
critical pressure. To verify this assertion, we ob- 
serve that the potential energy V is a functional 
of the displacement components (wv, v,w) and the 
pressure p. Let Class I be the set of all continuous 
differentiable functions (u,v, w) that satisfy the 
forced boundary conditions. 

It has been found that there exists a pressure 


“p’ for a given shell, such that a buckled configura- 


tion will persist if p >p’. The buckled configura- 
tion, being stable, provides a relative minimum 
to V among functions of Class I. This minimum 
of V depends only on p; hence, it will be denoted 
by f(p). (See Fig. 4.) 

It has been found that there exists a pressure 
p” (the Tsien critical pressure), such that f(p) > 0 
in the range p’ << p < p”, f(p”) = 0, and f(p) <0 
in the range p > p”. 

Let Class II be a given subset of Class I, as 
determined, for example, by assumptions about the 
nature of the deformation pattern. We have em- 
ployed two assumptions of this type: (1) the shell 
buckles without incremental hoop strain. (2) The 
function w, is represented by a single term of a 
Fourier series in # (see Eq. 36). When the func- 
tions (u,v, w) are restricted to Class II, the mini- 
mum value of V is ¢(p). If our assumptions are 
good, ¢(p) differs but slightly from f(p). 

Since Class IT is a subset of Class I, 6(p) 2f(p). 
Consequently, the graphs of f(p) and $(p) have 
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the general features shown by Fig. 10. The essen- 
tial characteristics of these functions are that ¢ 
and f are positive for small values of p and nega- 
tive for large values of p, and that the curve repre- 
senting ¢(p) lies above or in contact with the 
curve representing f(p). 

The Tsien critical pressure is the intercept of 
the graph of f(p) with the p-axis (Fig. 10). Evi- 
dently if ¢(p) is used as an approximation for f(p), 
the computed value of the Tsien critical pressure is 
too high. 


17. Potential Energy Barriers 


The maximum on the curve corresponding to p» 
(Fig. 4) represents a potential energy barrier that 
the shell must cross to arrive at the buckled form, 
if the pressure is exactly equal to the Tsien eriti- 
cal pressure. Therefore, it serves as a rough indica- 
tion of the imminence of snap-through. The value 
of this maximum may be derived from Eq. (100). 
For brevity, Eq. (100) is written as follows: 


y = ax — bx? + cz? 


AV se 
Oe Phy Ge 2 a 
Y EahL’? * Wo (a) 
a= B,— Kam, b= —(B.+ Ka), c= B3 


The graph of y is as shown in Fig. 11. 
The maximum value of y is determined by dif- 
ferentiation with respect to x. Thus, 
a — 2bx% + 3c2? = (b) 
The roots of Eq. (b) are 


ela oDe ay _ dae | 
eerie ek E : b? | 
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y The root x, provides the 
maximum, and the root v. 
x provides the minimum 
(Fig. 11). Since the value 


(¢) 


Fig. 11. Potential Energy 


Barrier of the minimum is zero, 
a — ba, + ca? = 0 (d) 
Eqs. (ce) and (d) yield 
b =2V/ ac (e) 
Consequently, Eq. (c) yields 
= TOES mvete Gio 
2 oh, ee NE (f) 


Eqs. (a) and (f) yield 


4 
Ymax = ae a9 (g) 
In terms of our previous notations, this equation 
yields, 
AV ) _ 4 AL 
Eee Pl OP 
where W, is the root of Eq. (114) that corresponds 
to the point of tangeney with the x-axis (Fig. 11). 
Eq. (115) is plotted in Fig. 13 for a/h = 1000, 
X = 7 and for n = 2 to 20. 


BW (115) 


18. Numerical Example 


Consider a shell with the following proportions: 
a/h = 100, L/a = 0.6010. The value L/a = 0.6010 
is selected to coincide with a tabulated value. This 
condition is unimportant. If a selected value of 
L/a does not appear in the tables, interpolation 
must be used. 


(a). Huler Critical Pressure for Shell with Simply- 
Supported Ends and No Axial End Constraint. 
The equilibrium pressure corresponding to any 
given state of deformation is represented in the fol- 
lowing form: p= KEh/a. The constant K is evi- 
dently equal to the compression hoop strain that 
exists in the unbuckled shell at pressure p. The 
value of K corresponding to the Euler critical pres- 
sure is denoted by AK;. By Eq. (113), K; = B,/a,, 
where B, = b, + c,h*/a?. Accordingly, in this ex- 
ample, B, = b, + 0.0001 c,. The constants a, bi, 
c, for a shell with simply supported ends have been 
tabulated (Tables 2 to 20). The number of lobes 
in the buckled form must be determined by trial 
to minimize K;. For very long shells, n= 2. In 
general, n increases with decreasing L or h. In the 
present example, L/a is small, but h/a is large. 
Therefore, a moderate value of n—for example, 
n = 10 — might be estimated. It is found by several 
trials that the value n=9 actually provides a 
minimum to K;. For n= 9, Table 9 yields (with 
AK=7) a, = 0.9327, b; = 0.0006329, c, = 10.45. 
Consequently, 8B, = 0.001678. Accordingly, Eq. 
(113) yields AK; = 0.001678/0.9327 = 0.001799. 
The condition X= 7 indicates that a uniform 
hydrostatic pressure acts on the ends of the shell. 
If \ = 0, the axial force due to the pressure on the 
ends is removed. Then a, = 0.7952, as noted at the 
bottom of Table 9. Since b, and c, are independent 
of A, B, has the same value as before. Thus, if 
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X = 0, then K; = B,/a, = 0.001678/0.7952 = 
0.002110. Accordingly, in this example, the hydro- 
static pressure on the ends reduces the Euler criti- 
eal pressure about 15%. 

For the case A = 0, von Mises derived a for- 
mula that may be put in the following form: 


1 he 
a pa EMO) SS ee 
he y(U ef: ae ) LZ ey) 


2n?—1—y 

ieee aA oT? 
De cas 

wa 


In the present numerical example, von Mises’ for- 
mula yields K; =0.001900. This result is about 10% 
lower than that computed by the present theory. 
There is seemingly a systematic deviation between 
the present infinitesimal theory and von Mises’ 
theory, for thick shells that are short compared to 
their radii. However, in all cases, the Tsien crite- 
rion yields lower values than the von Mises’ theory. 
(b) Presswre-Deflection Curves for Elastic Shell 
with Simply-Supported Ends and No Axial End 
Constraint. 


(116) 


The pressure-deflection curve is essentially a 
graph of K versus Wo, where K is defined as above. 
This curve may be plotted by means of Eq. (111). 
Setting n=9, we obtain from Table 9 (with A=7), 


a = 0.93827, a, = —1.270, 
b; = 0.0006329, b: = —0.02334, 63; = 0.3077, 
1 = 10.45, (Gy a 14.89, (AS 0.1124. 


The b’s and c’s are independent of ». The quanti- 
ties B,, Bz, B; are determined by B, = b, + ¢,h?/a’, 
Bz = bs +¢2 h?/a?, Bz = b3 + cz h?/a?. In the pres- 
ent example, h?/a? = 0.0001. Hence, 

B, = 0.001678, B. = — 0.02185, B, = 0.3077. 
Substituting these values of a1, d2, By, Bo, B; into 
Eq. (111), we obtain an equation whose graph is 
shown in Fig. 12. Since the curves corresponding to 
n = 8 and n = 9 intersect each other, the curve for 
n = 8 is also plotted. 

The intercept of the curve for n = 9 with the 
vertical axis is the Euler critical hoop strain, 
K;, = 0.001799. The minimum value on the curve 
for n= 9 is Kwin = 0.001170. The minimum value 
on the curve for n=8 is Kymin = 0.001080. This 
is the lowest minimum that occurs for any value of 
n. Therefore, it determines the lowest pressure at 
which the shell will maintain a buckled form, if it 
is perfectly elastic. 
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Fig. 12. Buckling Coefficient K versus 
Deflection Parameter W, 


Hence, in this example, the lowest pressure at 
which the elastic shell will maintain a buckled 
form is 60% of the Euler critical pressure. 

The value of K corresponding to the Tsien criti- 
cal pressure has been denoted by K,,. For this case, 
Ky, = 0.00133. This result may be obtained from 
Eq. (112). The Tsien eritical pressures have been 
marked on the curves of Fig. 15. The Euler critical 
pressure for this case (see above) is K; = 0.001799 
which is 35% higher than K,,;. 


(c) Effect of Rigid Ends. 

If the ends of the shell are hinged, but the axial 
displacements are constrained by the action of 
rigid end plates, the buckling pressure is increased 
significantly. Eqs. (111) and (113) remain valid, 
but the coefficients b,, bs, bs are changed. The con- 
stants @,, Go, C1, Co, Cz are not altered. 

The constants b,, bs, bs; have not been tabulated 
for a shell with rigid ends. Consequently, their 
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values must be computed by means of Eqs. (105), 
(107, (110), and (72). 

The constraint imposed by rigid end plates 
generally increases the number of lobes in the 
buckled form. Trying n= 10, we obtain by Eq. 
(48), €= 1.778. Hence, by Eq. (105), 

F,(é) = —0.003675, 
F3(€) = —0.002098. 


Using the values of the K’s from Table 1, we ob- 
tain by Eq. (107), 


F2(é) = 0.001648, 


vi = —0.0001565, Ye» = 0.005605, 
v3 = 0.000382, vs = —0.05900, 
Ws = 0.02203, vs = —0.006755. 


Eq. (110) yields X, = 0.002364. 

In Eq. (72), the functions fy, fo, .. . are to be re- 
placed by wn, ve, ..., respectively. Also, X, replaces 
¢;. The functions ¢. and g are discarded since they 
are negligible. 


Hence, 6, = 0.0006858, b. = — 0.01210, b; = © 
0.2929. Interpolating values from Table 10, we 
obtain 

a, = 0.904, a2 = — 1.27, c, = 11.68, c. = 15.14, 
cz; = 0.139. : 
Since B, = b, + ¢,h?/a?, B= be-+ 6,h?/a?, and 
Bz = bz + c3h?/a?, B, = 0.00185, B, = — 0.0106, 
B, = 0.293. Eq. (111) now provides a graph of K 
versus W, (Fig. 12). 

It is necessary to repeat the calculations for 
several other values of n. It is found by trial that 
the value n = 11 provides the lowest buckling pres- 
sure. The curves corresponding to n=10 and 
n = 11 are plotted in Fig. 12. It is seen that these 
curves are significantly higher than the curves 
obtained for a shell without axial end constraints. 
Similar calculations have been performed for 
L/a = 1.159 and the results have been plotted in 
Fig. 12. All the curves for L/a = 1.159 are lower 
than the corresponding curves for L/a = 0.6010. 


V. SUMMARY 


A theory, based on an energy analysis, has been 
developed for the snap-through and post-buckling 
behavior of simply-supported ideal shells under 
the action of external pressure. The principal re- 
sults of the theory are given: (a) by Eqs. (71), 
(72), (97), (100), (101), (111), (112), and (113) 
for elastic shells whose ends are free to warp out 
of their planes, and (b) by Eqs. (72), (105), (107), 
(110), (111), (112), and (113), and the modifica- 
tions indicated in Article 13 for elastic shells whose 
ends are rigid plates. 


The main results of the computations are pre- 
sented in the form of tables and graphs. Tables 
1 to 20 list the parameters needed for calculation 
of the buckling coefficient K given by p., = KEh/a. 
The use of Tables 1 to 20 is illustrated by a nu- 
merical example (Article 18). Table 21 gives values 
of K, for elastic shells whose ends are free to warp 
out of their planes, as determined by the infinites- 
imal theory and the Tsien snap-through theory for 
various values of L/a and X, and for a/h = 1000. 
For no axial pressure (A = 0), some values of K 


Legend | 


a/h=/00 
/nfinitesimal theory 


/ 

2 a/h=1000 
3 a/h=100 
4 


a/h=/000 } Elastic snap-through theory 


Log, K-10” 


=La | 


DNS 


Fig. 13. 


Log,/OL/a 


31 


Buckling Coefficients for Cylindrical Shells Subjected to Hydrostatic Pressure 


32 ILLINOIS ENGINEERING EXPERIMENT STATION 


| Oh =/000, A=7 


ae 


+ 


/0 25 ESO IO. ao BO HO 
L/a 


Fig. 14. Potential Energy Barriers Separating Buckled 
and Unbuckled Forms 


as calculated by von Mises’ theory are given for 
comparison. Table 22 lists similar values of K for 
a/h = 100. 

Discrepancies between von Mises’ theory and 
the present infinitesimal theory are greatest for 
short thick shells. Apparently, the trouble hes in 
the assumption that the shell buckles without in- 
cremental hoop strain. Von Mises did not make 
this assumption. 

For elastic cylinders whose ends are free to 
warp out of their planes, the Euler buckling coeffi- 
cient (Infinitesimal Theory) and the Tsien buck- 
ling coefficient (Elastic Snap-Through Theory) are 
plotted versus Z/a in Fig. 13 for a/h = 100, and 
a/h = 1000 with AX = z. Some of the data of Fig. 
13 are reproduced in Tables 21 and 22. For long 


slender cylinders (see Tables 21 and 22), the Euler 
buckling coefficient is only slightly higher than the 
Tsien coefficient. However, for relatively small 
values of L/a (say, L/a = 0.6010), the Euler coeffi- 
cient may be 30 to 35% higher than the Tsien 
coefficient. In the numerical example of Article 18, 
the Tsien coefficient is approximately 14% higher 
than the minimum pressure under which the elastic 
shell will maintain a buckled form. Prevention of 
end warping raises the critical pressure (Fig. 12). 

For 4=0 (no end pressure), all the critical 
pressures as determined by the different theories 
are raised. The effect of axial compression is great- 
est for small values of L/a. It becomes insignificant 
for very large values of L/a (Tables 21 and 22). 

The negative slopes of the load-deflection curves 
(Fig. 12) denote a condition favorable to snap- 
through. The potential energy barrier that the shell 
must overcome to snap-through is discussed in 
Article 17. Fig. 14 is a chart that shows these bar- 
riers for a/h = 1000 and A = z. The curve is dis- 
continuous because of sudden changes in n. The 
dashed curves have no significance; they merely 
outline the region in which the discontinuous curve 
lies. The points of discontinuity correspond to the 
cusps on curve 4 of Fig. 13. For example, if L/a = 
0.6 and a/h = 1000, Fig. 14 shows that n= 17 
and 10° A V/Ea* = 11.5. Hence, if a = 20 in. and 
E = 30,000,000 psi, AV = 2.76 in.-lb = 0.23 ft-lb. 
This result means that only 0.23 ft-lb of work must 
be supplied from the outside to cause snap-through. 
Accidental disturbances might easily supply this 
much energy. Imperfections are perhaps a more 
frequent cause of snap-through than accidental dis- 
turbances, although submarine hulls may be sub- 
jected to damaging shocks. 
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Ki 
.43787 
.061575 
-017515 
0268417 
0232170 
- 0217104 
0399290 
-0361575 
0340208 
- 0327367 
0319271 
0313963 
.0310364 
0178540 
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-0147512 
-0137773 
- 0130408 
0424754 
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Kio 
26572 
. 40644 
- 54541 
.68375 
.82178 
. 95964 
.09740 
.2351 
.38727 
.5103 
6479 
7854 
.9230 
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-7480 
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Table 1 


Values of K's for ¥ = 0.30 


Subscripts denote number of zeros preceding first significant figure. 


6 
If \=0, Kaui = 1.0472 and a2=0.07118. The b's and c's are independent of \. 
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Ko 
33126 
046584 
-013250 
.0:51760 
. 0024338 
- 0212940 
0375116 
0346584 
0530419 
0320704 
0314580 
-0310563 
- 0478405 
0559418 
0545849 
0135944 
0428576 
0123004 
0418727 
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Kis 
-027489 
-017671 
.013090 
010412 
- 0286504 
0274009 
0.64680 
0257446 
0251671 
- 0246953 
- 0243026 
0239706 
0236862 
- 0234399 
0232245 
0930845 
0228656 
0227146 
0225787 
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vy = 0,30 
bs 

039558 
034230 
032144 
051197 
017193, 
OWB041 
0.1494 
056222 
051996 
058235 
058989 
- 051270 
075228 
072535 
071046 
083360 
081408 
096960 
Owl98 
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Ks Ks Ks 
0.0: 1.44334 0.68539 0.16004 
0. 0.40657 0.19277 0. 
0. le 0.19723 0.093213 0. 
0. 0.26701 0.11807 0.055688 0. 
0.012798 0.079100 0.037263 0.1 
0 0568822 0.056878 0.026773 0. 
0. 0540245 0.042946 0.020204 0. 
0 0525084 0.033611 0.015807 0.15445 
0.0516438 0.027041 0.012713 0.15441 
0.0511218 0.022236 0.010452 0.15437 
0.0679152 0.018613 0. 0287479 0.15435 
0. 0657437 0.015813 0. 0274309 0.15433 
0.0642685 0.013603 0.0: woes 0.15431 
0 .0632380 0.011827 0.025556 0.15430 
0 .0625006 0.010379 0. 0.15429 
0 .0619617 0. Oug 0.02 0, 15428 
0.0515605 0. 0. 0.15427 
0.0612568 Oe 0. 0.15427 
0.0610235 0. 0.0:31075 0.15426 
Table 1 (Concluded) 
Kis Kis Ku 
0.43982 2.5839 0.086128 
0.33576 0.09895 0.020763 
0.30718 0.51677 0. 0282027 
0.29502 0.32298 0 0240780 
0.28868 0.22147 0.0223231 
. 16035 0.28493 0.16149 0.0.14493 
13963 0.28254 0.12304 0 ..0396502 
12370 0.28091 0.096894 0. 0367496 
.11107 0.27975 0.078299 0. 0349059 
10079 0.27890 0.064596 0 0336778 
.092271 0.27826 0.054207 0 .0328282 
- 085085 0.27776 0.046140 0 .0322216 
.078943 0.27736 0.039752 0.0317769 
.0736381 0.27704 0.034605 00314435 
068992 0.27678 0.030398 0..0311886 
.064905 0.27656 0.026915 0.0199038 
.061276 0.27638 0.023999 0..0183393 
.058032 0.27623 0.021532 00370878 
.055116 0.27610 0.019427 0 0460749 
Table 2 
Values of Coefficients for Computing Buckling Loads —n = 2 
n=2 =r 
a az bi bs 
1.0943 0.06018 0.039680 
1.0773 0.06414 0.034165 —0. 037060 
1.0681 0.06629 0.082065 (i) 50: 
1.0625 0.06759 0.031134 —0.05 
1.0589 0.06843 0.046720 —0.0; 
1.0547 0.06942 0.042788 —0.05 
1.0524 0.06996 0.041854 =. 
1.0506 0.07040 0.055574 —=()), 
1.0491 0.07074 0.051769 —0. 05° 
1.0484 0.07090 0. 0672388 —0.051254 
1.0480 0.07098 0 098483 —0.0,6066 
1.0477 0.07107 0.051094 —0. 061921 
1.0475 0.07111 0.074439 —0.0;7986 
1.0474 0.07113 0.072111 —0.073838 
1.0473 0.07115 0.038454. —0.071617 
1.0473 0.07116 0.052501 —0. 085054 
1.0472 0.07117 0.009686 —0.032122 
1.0472 0.07117 0.004119 —0.0s1011 
1.0472 0.07117 0.009321 — 0.093034 
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Ky 

. 21035 
-084115 
.045660 
028753 
.O19794 
.014467 
.011039 
0287018 
. 0270368 
- 0258084 
0248762 
0241519 
-0235779 
0231153 
0227366 
0224237 
0221613 
0219394 
0217500 


Kis 


.04720 
88357 
. 83776 
-81812 
.80784 
-80176 
79786 
- 79522 
79333 
-79194 
.79089 
. 79007 
.78943 
.78890 
.78848 
- 78812 


78783 


-78758 
78736 


cr 

3233 
.3102 
8032 
-2990 
- 2963 
- 2932 
2915 
2901 
2890 
. 2886 
2883 
2880 
2879 
2878 
. 2878 
2877 
2877 
2877 
2877 
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Ks 
087164 
090544 
.091700 
. 092232 
.092518 
.092691 
. 092804 
092880 
-092935 
.092975 
093006 
093030 
- 093049 
.093064 
.093076 
093087 
. 093096 
.093103 
.093110 
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Ku 


071177 
-072579 
.073047 
.073259 
073374 
073442 
- 073487 
073517 
073539 
.073555 
073567 
073576 
.073584 
.073590 
-073595 
073599 
.073603 
.073606 
-073608 


C2 
05798 
05634 
05553 
-05508 
05480 
05448 
05431 
-O5417 
05407 
05402 
05400 
05397 
.05396 
.05396 
05395 
05395 
053894 
05394 
.05394 
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Ks 
-8422 
-5157 
-0582 
-81969 
.67139 
56956 
49506 
43805 
39296 
35638 
.32609 
30058 
. 27879 
. 25997 
- 24354 
. 22908 
.21624 
20477 
.19446 
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Kis 


0.25839 
0.11074 
0.061520 
0.039149 
0.027103 
0.019876 
0.015199 
0.011999 
00297137 
0 ..0280244 
0.0267405 
0 0257419 
0 0249499 
0.243112 
0 .0237886 
00233556 
0..0229929 
0..0226859 
00224239 


cs 


- 025626 
. 025624 
025622 
025621 
025621 
025620 
025620 
025620 
025619 
-025619 
025619 
025619 
025619 
- 025619 
025619 
025619 
- 025619 
025619 
-025619 
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Table 3 
| Values of Coefficients for Computing Buckling Loads —n = 3 
: n=3 =r v = 0,30 
FS L/a a a2 bi be bs C1 c2 3 
; 2.4 2.704 0.9664 0.025343 0.021910 —0.027619 0.01077 0.8671 0.1730 0.01262 
2.8 3.155 0.9444 0.02318 0.021120 —0.024485 0. 026490 0.8052 0.1519 0.01262 
3.2 3.606 0.9301 0.03476 0.036956 —0.022793 0.024137 0. 0.1424 0.01263 
4 4.508 0.9134 0.04837 0.033063 —0.021234 0.021911 0. 0.13824 0.01263 
5 5.634 0.9026 0.05709 0.031318 — 0.035394 0.038630 0. 0.1275 0.01264 
6 6.761 0.8968 0.06182 0.016536 —0.032639 0 .0s4443 0.6 0.1252 0.01264 
th 7.888 0.8933 0.06468 0.043589 —0.031450 0. 0. 0.1240 0.01264 
8 9.015 0.8910 0.06653 0.042128 — 0.018596 0.03152 0. 0.1233 0.01264 
10 11.269 0.8883 0.06871 0.058834 —0.043569 0. 0:68 506 0.68 0.1226 0.01264 
12 13.523 0.8869 0.06989 0.054292 —0.011754 0.04322 0. 0.1222 0.01264 
15 16.904 0.8857 0.07086 0.051768 —0.057145 0.0: 0. 0.1219 0.01264 
20 22.538 0.8848 0.07161 0.095622 —0.052271 0.054352 0.6 0.1216 0.01264 
25 28.172 0.8843 0.07196 0.052307 —0.069335 0.051799 0. 0.1216 0.01264 
30 33.81 0.8841 0.07215 0.061113 —0. 064494 0.068725 0. 0.1215 0.01264 
40 45.08 0.8839 0.07234 0.073527 —0.061421 0.062780 0. 0.1214 0.01264 
50 56.34 0.8838 0.07242 0.071454 —0.07583 par eeae 0. 0.1214 0.01264 
60 67.61 0.8837 0.07247 0.036938 —0.07279 0. 0.6 0.1214 0.01264 
75 84.52 0.8836 0.07251 0.082879 07 0. 0.6 0.1214 0.01264 
100 112.69 0.8836 0.07254 0.099986 —0.0534 0. 0.6 0.1214 0.01264 
If A\=0, a1=0.8836 and az=0.07258. The b’s and c’s are independent of A. 
Table 4 
Values of Coefficients for Computing Buckling Loads —n = 4 
n=4 N=9r pv = 0.30 
fa L/a a az bi be bs al c2 
2.0 1.690 0.9402 —0.09791 0.021769 —0.01273 0.03261 1.719 0.4779 0. 
2.4 2.028 0.9089 —0.04568 0.039656 —0.026993 0.01837 1.532 0.3522 0. 
2.8 2.366 0.8900 —0.01418 0. —0 004116 0.01112 1.425 0.2959 0. 
3.2 2.704 0.8778 0.026266 0. —0 . 022562 0.027116 1.358 0.2671 0. 
4 3.381 0.8634 0.03031 0.051548 —0.021132 0.023311 1.281 0.2410 0.02251 
5 4.226 0.8542 0.04569 0.016660 —0.034877 0.021507 1,233 0.2289 0.02250 
6 5.071 0.8492 0.05405 0.043302 — 0.032421 0.037797 1.208 0, 2237 0.02249 
i 5.916 0.8461 0.05909 0.041812 —0.031330 0.034422 1.192 0.2211 0.02248 
8 6.761 0.8442 0.06236 0.041074 —0.017886 0.0; 1.183 0.2196 0.02248 
10 8.452 0.8419 0.06621 0.054454 —0.013275 0.03115 u ei by A § 0.2180 0.0% 
12 10.14 0.8406 0.06830 0.052161 —0.041592 0.045 1.165 0.2172 0. 
15 12.68 0.8396 0.07001 0, 068881 —0.05;6562 0 1.160 0.2167 OF 
20 16.90 0.8388 0.07134 0.052809 —0.052088 0 1.156 0.2163 0.02248 
25 21.13 0.8384 0.07195 0.051145 —0.068585 0 1.154 0.2161 0.02248 
30 25.35 0.8382 0.07229 0.075478 —0.064152 0 1.153 0.2160 0.02248 
40 neceh 0.8380 Beane 0.071698 —0.061322 0 1.152 0.2159 0.02248 
0 42.2 0.8379 0.07277 0.036805 —0.07545 0 1.152 0.2158 0.02248 
If A\=0, ee 8378 and a2=0.07305. The b's and c’s are independent of 2. 
Table 5 
Values of Coefficients for Computing Buckling Loads —n— 5 
n=5 \=r v = 0.30 
é L/a a az br be bs C1 cz C3 
1.6 1.082 0.9676 —0.3429 0.02170 —0.02441 0.09734 3.248 1.731 0.03476 
2.0 1.352 0.9138 —0.1931 0.021080 —0.01227 0.04999 2.675 0.9347 0.03489 
2.4 1.623 0.8846 —0.1117 0.035895 —0.026733 0.02822 2.388 0.6409 0. 
2.8 1.893 0.8669 —0.06262 0.033457 —0.023963 0.01712 2,223 0.5111 0. 
3.2 2.164 0.8555 —0.03077 0.052147 —0.022467 0.01096 2.119 0.4459 0. 
4 2.704 0.8420 0.026678 0.019454 —0.021089 0.025124 2.000 0.3883 0 
5 3.381 0.8334 0.03065 0.014068 —0.034693 0. : 1.926 0.3623 0. 
6 4.057 0.8288 0.04367 0.042017 —0.032329 0. 4 1.886 0.3518 0. 
7 4.733 0.8259 0.05151 0.011108 —0.031280 0.036890 1.863 0.3466 0. 
8 5.409 0.8241 0.05661 0..0;6567 —0.017586 0.034191 1.848 0.3438 0. 
10 6.761 0.8219 0.06261 0.052726 —0.013140 0.051802 1.830 0.3409 0. 
12 8.113 0.8208 0.06686 0.051324 —0.041530 0. 048942 1.820 0.3396 0. 
15 10.142 0.8198 0.06852 0.065456 —0.0;6305 0.053757 1.812 0.3386 0. 
20 13.52 0.8191 0.07060 0.061734 —0.052005 0.041215 1.806 0.3379 0. 
25 16.90 0.8187 0.07155 0.077115 —0.068229 0.055027 1,803 0.3376 0. 
30 20.28 0.8185 0.07208 0.073432 —0.064012 0.052439 1.802 0.3375 0. 
40 27.04 0.8184 0.07259 0.071086 —0.051257 0.067783 1.800 0.3373 0. 
50 33,81 0.8183 0.07283 0.084452 —0.075173 0.053205 1.799 0.3372 0. 
If \=0, ai=0.8181 and a2:=0.07326. The b's and c's are independent of d. 
Table 6 
Values of Coefficients for Computing Buckling Loads —n— 6 
n=6 N= vy = 0.30 
— L/a a az by be bs C1 C2 C3 
1.2 0.6761 1.0639 —0.9902 0.023362 —0.05376 0.3080 6.707 9.160 0.04956 
1.6 0.9015 0.9519 —0.5249 0.021470 —0.02393 0.1386 4.663 3.303 0.05000 
2.0 1.127 0.9000 —0.3095 0.087312 —0.01202 0.07128 3,844 1.677 0.05021 
2.4 1.352 0.8719 —0.1925 0.033991 —0 026542 0.04028 3.434 1.083 0.05032 
2.8 1.578 0.8549 —0.1220 0.032341 —0.023888 0.02445 3.198 0.8222 0.05038 
3.2 1.803 0.8438 —0.07619 0.031454 —0. 022418 0.01570 3.049 0.6925 0.05043 
4 2.254 0.8309 —0.02234 0.046401 —0.021067 0.007341 2.879 0.5798 0. 
5 2.817 0.8226 0.01111 0.042754 —0.034600 0.023356 2.773 0.5302 0. 
6 3.381 0.8181 0.03083 0.041365 —0.032282 0.001743 2.716 0.5106 0. 
7 3.944 0.8154 0.04212 0.057497 —0.031254 0.039908 2.682 0.5014 0. 
8 4.507 0.8136 0.04944 0.054444 —0.0:7385 0.036030 2.660 0.4963 0. 
10 5.634 0.8115 0.05806 0.051844 —0.0s3087 0.032594 2.634 0.4914 0. 
12 6.761 0.8104 0.06274 0.058956 —0.031500 0.031288 2.621 0.4892 0. 
15 8.452 0.8095 0.06657 0.053687 —0.056181 0.015412 2.609 0.4877 0. 
20 11.27 0.8088 0.06954 0.061170 —0.051966 0.011749 2,600 0.4866 0. 
25 14.08 0, 8084 0.07092 0.074790 —0.068035 0.057237 2.596 0.4862 0. 


If \=0, a1=0.8078 and a:=0.07337. The b’s and c’s are independent of d. 
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Table 7 
Values of Coefficients for Computing Buckling Loads —n = 7 
n=7 A=c7 v = 0.30 
é L/a ay az bi be bs 
1.0 0.4829 1.1624 —2.009 0.023948 —0.08504 0.6668 
1.2 0.5795 1.0522 —1.373 0.022433 —0.05314 0.4162 
1.6 0.7727 0.9426 —0.7400 0.021064 —0.02364 0.1875 
2.0 0.9659 0.8919 —0.4472 0.035292 —0.01188 0.09646 
2.4 1.159 0.8644 —0.2881 0.032889 —0. 026580 0.05455 
2.8 1.352 0.8478 —0.1922 0.031694 — 0.023838 0.03314 
3.2 1.545 0.8370 —0.1299 0.031052 — 0.022388 0.02129 
4 1.932 0.8243 —0.05671 0.014634 —0.021054 0.029965 
5 2.415 0.8162 —0. 029858 0.051994 —0. 034544 0.024560 
6 2.898 0.8118 0.01559 0.0;9889 —0.032255 0.022370 
7 3.381 0.8091 0.03094 0.055431 —0.031239 0.021348 
8 3.864 0.8074 0.04090 Mane ECs — 0.037345 0.038206 
10 4.829 0.8054 0.05262 0.051330 —0.0s3050 0.033532 
If \=0, a1=0.8018 and az=0.07344. The b's and c’s are independent of 2. 
Table 8 
Values of Coefficients for Computing Buckling Loads — n= 8 
n=8 A=7 vy = 0.30 
— L/a a a2 by bz bs 
1.0 0.4226 1.153 —2.645 0.022993 —0.08441 0.8667 
1.2 0.5071 1.045 —1.814 0.001844 —0.05273 0.5411 
1.6 0.6761 0.9367 —0.9884 0.038064 —0.02346 0.2438 
2.0 0.8867 —0.6061 0.034012 —0.01178 0.1255 
2.4 0.8596 —0.3985 0.032190 —0.026470 0.07103 
2.8 0.8432 —0.2732 0.031284 — 0.023808 0.04316 
3.2 0.8326 —0.1920 0.017975 — 0.022370 0.02774 
4 0.8201 —0.09641 0.033512 —0.021046 0.01299 
5 0.8121 —0.03525 0.041511 —0. 034509 0.025950 
6 0.8077 —0.022024 0.057492 — 0.032237 0.023093 
4 0.8051 0.01801 0.054113 —0.0531213 0.021760 
8 0.8034 ee 0.052438 —0.0s7285 0.021072 
10 4. "296 0.8014 0.04630 0.051012 —0.013024 0.034614 
If A=0, a1=0.7979 and az=0.07349. The b’s and c’s are independent of 2. 
Table 9 
Values of Coefficients for Computing Buckling Loads —n = 9 
n=9 = v = 0.30 
E a ae bi bs bs 
1.0 a: 1.147 —3.366 0.002349 —0.08398 1.0934 
1.2 , 1.040 —=2 315 0.021448 —0.05292 0.6827 
1.6 6 0.9327 —1.270 0.036329 —0.02334 0.3077 
2.0 0. 0.8832 —0.7862 0.058149 —0.01172 0.1585 
2.4 0.¢ 0.8563 —0.5235 0.031719 —0 026437 0.08970 
2.8 i Me 0.8401 —0.3651 0.031008 —0.023788 0.05453 
3.2 step 0.8296 —0.2623 0.016260 —0. 022357 0.03505 
4 1.é 0.8172 —0.1414 0.042756 —0.021041 0.01642 
5 ie 0.8093 —0.06404 0.011186 —0.034484 0.027525 
6 2.3 0.8050 —0.02201 0.055880 —0.032225 0.023913 
7 2. 0.8024 0.023333 0.053228 —0.031222 0.022223 
8 3. 0.8007 0.01978 0.051914 —0.0s7247 0.021356 
10 3. 0.7987 0.03913 0.067941 —0.043008 0.035839 
12 4.6 0.7977 0.04964 0.058856 —0.0)1462 0.082901 
15 5. 0.7968 0.05823 0.091587 —0.0;6022 0.031068 
20 7.512 0.7961 0.06492 0.075037 —0.051914 0.0s3941 
If A=0, ai =0.7952 and az=0.07352. The b's and c’s are independent of A. 
Table 10 
Values of Coefficients for Computing Buckling Loads — n= 10 
= 10 =. vy = 0.30 
& L/a a ag bi bs bs 
0.8 0.2704 1.339 —6.559 0.023302 —0.1481 2.3217 
1.0 0.3381 1.143 —4.171 0.021894 —0.08367 1.3467 
1.2 0.4057 1.036 —2.874 0.021167 —0.05227 0.8410 
1.6 0.5409 0.9298 —1.584 0.035103 —0.02326 0.3792 
2.0 0.6761 0. 8807 —0.9876 0.032539 —0.01168 0.1953 
2.4 0.8113 0.8540 —0.6634 0.051386 —0.026412 0.1106 
2.8 0.9466 0.8379 —0.4677 0.048127 —0.023794 0.06724 
3.2 1.082 0.8275 —0.3410 0.0;5047 —0.022348 0.04323 
4 1.352 0.8152 —0.1918 0.042222 —0.021037 0.02026 
5 1.690 0.8073 —0.09624 0.059563 —0.0:4468 0.029287 
6 2.028 0.8030 —0.04437 0.094742 —0.032217 0.04831 
7 2.366 0.8005 —0.01309 0.052604 —0.031218 0.002745 
8 2.704 0.7988 0.027217 0.051543 —0. 017221 0.001675 
10 3.381 0.7968 0.03109 0.006407 —0.053000 0.087215 
12 4.057 0.7958 0.04406 0.063112 —0.0s1454 0.053586 
15 5. Po 0.7949 nee 5467 0.061282 —0.0;6006 0.031508 
20 6.7 0.7942 0.06293 0.0;4074 —0.051911 0.0:4882 
If \=0, a= if 7933 and az=0.07354. The b's and c’s are independent of 2. 
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Table 11 
Values of Coefficients for Computing Buckling Loads —n = 11 
=11 A= = 0.30 
— L/a a a2 by be bs 
0.7 0.2151 1.501 —10.40 0 023732 —0.1937 3.830 
0.8 0.2459 1.335 —7.949 0.002720 —0.1428 2.805 
1.0 0.3073 1.140 —5.061 0.021560 —0.08234 1.627 
1.2 0.3688 1.033 —3.492 0.039612 —0.05213 1.016 
1.6 0.4917 0.9277 —1.932 0.034202 — 0.02319 0.4581 
2.0 0.6147 0.8788 =1.210 0.032091 —0.01165 0.2360 
2.4 0.7376 0.8523 —0.8179 0.051141 —0.026394 0.1336 
2.8 0.8605 0.8363 —0.5814 0.016929 — 0.023763 0.08127 
3.2 0.9834 0.8259 —0.4279 0.014156 — 0.002342 0.05226 
4 1.229 0.8137 —0.2474 0.041830 —0.021034 0.02450 
5 1.537 0.8058 —0.1318 0.057876 — 0.034456 0.01123 
6 1.844 0.8016 —0.06907 0.053906 —0.032211 0.025844 
76 2.151 0.7990 —0.038123 0.052145 —0.031215 0.023327 
8 2.459 0.7974 —0.01667 0.051271 —0. 037204 0.022027 
10 3.073 0.7954 0.02221 0.05279 =e anki 0.038731 
12 3.688 0.7944 0.03790 0.062565 —0.0.1 0.034340 
15 4.610 0.7935 0.05074 0.061057 =(0)3 0.5004 0.031825 
20 6.147 0.7928 0.06072 0.073361 —0.051907 0.045906 
If \=0, a1=0.7919 and az=0.07356. The b’s and c’s are independent of \. 
Table 12 
Values of Coefficients for Computing Buckling Loads — n = 
n= 12 A\=7 vy = 0.30 
é L/a a ay br be bs 
0.7 0.1972 1.4976 —12.395 0.023128 —0.1933 4,552 
0.8 0.2254 1.3320 —9.4730 0.022279 —0.1425 3.333 
1.0 0.2817 1.1372 —6.0362 0.021307 —0.08328 1.933 
1.2 0.3381 1.0314 —4.1693 0.038056 —0.05204 1.208 
1.6 0.4508 0.9262 —2.3131 0.033522 —0.02314 0.5446 
2.0 0.5634 0.8775 —1.4539 0.031752 —0.01162 0. 2806 
2.4 0.6761 0.8510 —0.9872 0.019566 —0.026381 0.1589 
2.8 0.7888 0.8350 —0.7057 0.035610 —0 003756 0.09665 
3.2 0.9015 0.8247 —0.5231 0.013484 —0.022336 0.06216 
4.0 1.1269 0.8125 —0.3083 0.041534 —0.021032 0.02915 
5.0 1.4086 0.8047 —0.1708 0.036603 —0.034446 0.01336 
15.0 1.690 0.8005 —0.09615 0.053275 — 0.032207 On 0sg 258 
a0 1.972 0.7980 —0.05112 0.051799 —0.031212 0. 
8.0 2.254 0.7963 —0.02190 0.051066 —0.0:7187 0.0: 
Tf X=0, a1=0.79089 and avs=0.073567. The b's and c’s are independent of \. 
Table “13 
Values of Coefficients for Computing Buckling Loads — n= 13 
n= 13 AN=7 = 0.30 
£/ L/a a a2 bi be bs 
0.7 0.1820 1.4946 —14.558 0.022660 —0.1934 5.336 
0.8 0. 2080 1.3295 —11.1289 0.021938 —0.1423 3.907 
1.0 0.2600 1.1353 —7.0960 0.041112 —0.08443 2.267 
1.2 0.3121 1.0298 —4,9053 0.036850 —0.05195 1.416 
1.6 0.4161 0.9249 —2.7270 0.032994 —0.02310 0.6386 
2.0 0.5201 0.8764 —1.7188 0.031490 —0.01160 0.3291 
2.4 0.6241 0.8500 i se Ala 0 018132 —0.026381 0.1864 
2.8 0.7281 0.8341 —0.8409 0.054769 —0.023749 0.11384 
3.2 0.8322 0.8238 —0.6266 0.042962 —0 022333 0.07292 
4.0 1.040 0.8116 —0.3745 0.031304 —0.021030 0.03420 
5.0 1.300 0.8039 —0.2132 0.055611 —0.0;4439 0.01568 
6.0 1.560 0.7996 —0.1256 0.052782 —0.032203 0.028162 
7.0 1.820 0.7971 —0.07274 0.051528 —0.031210 0.04648 
8.0 2.080 0.7955 —0.03845 0.059055 —0.047177 0.022832 
If \=0, a1=0.79007 and a2=0.073576. The b’s and c’s are independent of , 
Table 14 
Values of Coefficients for Computing Buckling Loads —n — 14 
n= if N= vy = 0,30 
& L/a a a bi bo 
0.7 0.1690 1.4922 —16.8938 0.022289 —0.19283 6. 
0.8 0.1932 1.3275 —12.9170 0.021668 —0.14213 4: 
1.0 0.2415 1.1338 —8.2404 0.039568 —0.08305 2. 
1.2 0.2898 1.0286 —5.7000 0.035896 —0.05188 its 
1.6 0.3864 0.9239 —3.1741 0.032578 —0.02307 0. 
2:0 0.4829 0.8755 —2.0049 0.031282 —0.01159 0. 
2.4 0.5795 0.8492 —1.3698 0.017000 —0 .026363 0. 
2.8 0.6761 0.8334 —0.9869 0.054105 — 0.023744 0. 
3.2 0.7727 0.8230 —0.7383 0.042549 —0. 022330 0. 
4.0 0.9659 0.8110 —0.4460 0.041122 —0.021029 0. 
5.0 1.207 0.8032 —0.2590 0.054829 —0.034432 0. 
6.0 1.449 0.7990 —0.1574 0.052394 —0.032218 0. 
7.0 1.690 0.7965 —0.09609 0.051314 —0.031208 0. 
8.0 1.932 0.7948 —0.05632 0. 067788 —0.0s7162 0.0: 


Bul. 443. 


If \=0, a =0.78943 and a2=0,073584. The b’s and c’s are independent of A, 
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.2730 
2737 
- 2742 
2744 
2748 
2750 
.2751 
.2751 
2752 
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ILLINOIS ENGINEERING EXPERIMENT STATION 


If \=0, a1 =0.7878 and a2=0,07360. The b's and c's are independent of X. 


Table 15 
Values of Coefficients for Computing Buckling Loads —n = 15 
= 15 =7 vy = 0.30 
L/a ay a bi be bs C1 
0.1578 1.49031 —19.4030 0.0219918 —0.1926 7.0940 124.016 
0.1803 1.32592 — 14.8381 0. 0214514 —0.1420 5.1947 89.6418 
0.2254 1.13259 —9.4699 0 0383243 —0.08296 3.0135 56.3471 
0.2704 1.02757 —6.5538 0 .0351298 —0.05183 1.8823 41.4908 
0.3606 0.92315 —3.6543 0 .0322427 —0.02305 0.8492 28.9734 
0.4507 0.87482 —2.3123 0.0311157 —0.01158 0.4376 23 .9406 
0.5409 0.84857 — 1.5833 0 0460908 — 0.026356 0.2479 21.4086 
0.6310 0.83274 —1.1487 0.0435719 —0. 023740 0.1508 19.9506 
0.7212 0.82246 —0.85839 0. 0422182 —0.022328 0.09700 19.0320 
0.9015 0.81038 —0.52288 0.059769 —0.021028 0.04551 17.9794 
L127 0.80265 —0.30813 0.054204 — 0.034428 0.02087 17.3216 
1.352 0.79845 —0.19151 0. 0520846 —0.032197 0.01086 16.9695 
1.578 0.79591 —0.12118 0.0511447 —0.031207 0.026185 16.7589 
1.803 0.79427 —0.075527 0. 0867869 —0.037156 0.023768 16.6229 
2.254 0.79234 —0.021845 0.0528178 —0.0s2971 0.021624 16.4637 
2.704 0.79129 +0 .027316 0.0513692 —0.051443 0.038071 16.3776 
If \=0, a1=0.7889 and a2=0.07359. The b's and c’s are independent of 2. 
Table 16 
Values of Coefficients for Computing Buckling Loads —n = 16 
= 16 }=7 vy = 0.30 
L/a ay a2 bi be bs C1 
0.1479 1.4887 — 22.0856 0.021749 —0.1924 8.067 141.0 
0.1690 1.3246 — 16.8920 0.021274 —0.1419 5.907 101.9 
0.2113 1.1316 —10.7844 0.037308 —0.08288 3.427 64.09 
0.2535 1.0268 —7.4667 0.034503 —0.05178 2.141 47.20 
0.3381 0.9225 —4.1678 0.031969 —0.02303 0.9658 32.96 
0.4226 0.8743 —2.6409 0.039794 —0.01156 0.4978 27.24 
0.5071 0.8480 —1.8115 0.035346 —0.06349 0.2820 24.36 
0.5916 0.8322 —1.3114 0.043135 —0.093754 0.1715 22.70 
0.6761 0.8220 —0.9868 0.031947 —0.022340 0.1103 21.65 
0.8452 0.8099 —0.6050 0. 058574 —0.021035 0.05177 20.46 
1.056 0.8022 —0.3607 0.053689 —0.034466 0.02374 19.71 
1.268 0.7980 —0.2280 0.051829 —0.032194 0.01236 19.31 
1.479 0.7955 —0.1480 0.051004 —0.031205 0.027037 19.07 
1.690 0.7938 =0% 09606 0.055951 —0.037140 0.024288 18.94 
If \=0, a1=0.7885 and a2=0.07360. The b's and c’s are independent of 2. 
Table 17 
Values of Coefficients for Computing Buckling Loads — n= 17 
=17 =r vy = 0.30 
L/a a a2 bi bo bs cy 
0.1392 1.4874 — 24.9387 0.021548 —0.1923 9.100 159.14 
0.1591 1.3235 —19.0764 0.021128 —0.1418 6.664 115.05 
0.1989 1.1308 —12.1824 0.036468 —0.08283 3.867 72.34 
0.2386 1.0261 —8.4875 0.033986 —0.05175 2.416 53.27 
0.3182 0.9220 —4.7139 0.031742 —0.02302 1.090 37.21 
0.3977 0.8738 —2.9904 0.038669 —0.01156 0.5618 30.75 
0.4773 0.8476 —2.0542 0 .0s4732 —0.026347 0.3182 27.50 
0.5568 0.8318 —1.4897 0.052775 —0. 023734 0.1936 25.62 
0.6364 0.8216 —1.1233 0.011723 —0.022324 0.1245 24.44 
0.7954 0.8095 —0.6924 0.057589 —0.021026 0.05843 23.09 
0.9943 0.8018 —0.4166 0.053265 —0.034422 0.02680 22.25 
1.193 0.7976 —0.2668 0.051619 —0.032194 0.01395 21.80 
1.392 0. —0.1765 0. 068890 —0.031205 0.027944 21.52 
1.591 0.7 —0.1179 0.065269 —0.0:7148 0.0:4840 21.35 
If \=0, a1=0.7881 and a0 07360. The b's and c’s are independent of i. 
Table 18 
Values of Coefficients for Computing Buckling Loads — n= 18 
=18 =r vy = 0,30 
L/a a a2 bi bz bs C1 
0.1315 1.4864 — 27.967 0.021379 —0.1922 10.199 178.35 
0.1502 1.3227 —21.395 0.021005 —0.1417 7.469 128,94 
0.1878 1.1301 —13.666 0.035765 —0.08279 4.334 81.08 
0.2254 1.025 —9.468 0.053552 —0.05171 2.707 59.71 
0.3005 0: —5.294 0.031553 —0.02300 1.222 41.71 
0.3756 0. —3.361 0.037727 —0.01155 0.6296 34.47 
0.4507 0. ‘ —2.312 0.04218 —0.026342 0.3567 30.82 
0.5259 0.8315 —1.679 0.042473 —0.023732 0.2170 28.73 
0.6010 0, 8212 —1.2682 0.051536 —0.022325 0.1396 27.40 
0.7512 0.8092 —0.7851 0. 056764 —0.021026 0.06551 25.89 
0.9391 0.8015 —0.4760 0.052910 —0.034417 0.03005 24.94 
1.127 0.7973 —0.3081 0.051443 —0.032192 0.01564 24.44 
1.315 0.7948 —0.2068 0.067923 —0.031204 0. 0e889 1 24.13 
1.502 0.7932 —0.1411 0. 0¢4089 —0.047137 0.025428 23.94 


c2 


2835 .839 
1664 .536 
684.529 
332,244 
107 .647 
46.1196 
23.9515 
14.4137 
- 76469 
- 85008 
. 22540 
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- 24099 
12838 
-08521 
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C3 


294380 
299463 
305441 
308688 
311916 
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-315961 
315994 


C3 


3349 
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Table 19 
Values of Coefficients for Computing Buckling Loads —n = 19 
n=19 A}=r v = 0.30 
L/a a ay br be bs a ce c3 
. 1246 1.4854 —31.1693 0.021237 —0.1921 11.360 198.65 7292.1 0.4719 
. 1423 1.3219 — 23.8467 0.039015 —0.1416 8.319 143.63 4278.0 0.4802 
-1779 1.1295 —15.235 0.035171 —0.08275 4.827 90.32 1756.7 0.4899 
. 2135 1.0250 —10.5576 0.033186 —0.05169 3,016 66.52 850.6 0.4951 
. 2847 0.9212 —5.9065 0.031393 —0.02299 1.361 46.47 273.2 0.5004 
.3559 0.8731 —4.1669 0.016930 —0.01155 0.7014 38.40 115.0 0.5028 
4270 0.8470 —3.7536 0.033783 —0.0:6339 0.3973 34.34 58.27 0.5041 
-4982 0.8312 —2.5842 0.042219 —0.023731 0.2417 32.00 33.84 0.5049 
. 5694 0.8210 —1.4214 0.031378 —0. 022321 0.1555 30.53 21.94 0.5054 
or bilrd 0.8087 —0.8832 0.056067 —0.021025 0.07298 28.84 11.96 0.5060 
. 8896 0.8010 —0.5388 0.052611 —0.034417 0.03348 27.79 7.833 0.5064 
.068 0.7968 —0.3516 0.051294 —0.032192 0.01742 27.22 6.334 0.5066 
246 0.7943 —0.2388 0.067109 —0.051204 0.009905 26.89 5.682 0.5067 
-423 0.7927 —0.1556 0.054214 —0.0s7141 0.026047 26.67 5.361 0.5068 
0.7876 and az=0.07361. The b's and c’s are independent of ). 
Table 20 
Values of Coefficients for Computing Buckling Loads —n = 20 
n = 20 =r v = 0.30 
L/a a a2 bi be bs C1 cz cs 
1183 1.4847 —34.5428 0.021116 —0.1920 12.584 220.06 8950.8 0.52288 
. 1352 1.3213 — 26.4295 0.038132 —0.1416 9.215 195.11 5250.7 0.53203 
. 1690 1.1291 —16.8884 0.034664 —0.08272 5.347 100.06 2155.6 0.54278 
-2028 1.02464 —11.7056 0.032874 —0.05167 3.340 73.70 1043.3 0.54862 
-2704 0.9208 —6.5522 0.031256 —0.02298 1.507 51.48 334.6 0.55443 
.3381 0.8728 —4.1669 0.036251 —0.01154 0.7770 42.55 140.6 0.55712 
-4057 0.8467 —2.8712 0.043412 —0.026336 0.4402 38.05 70.87 0.55858 
.4733 0.8309 —2.0899 0.042001 —0.023692 0.2678 35.46 40.90 0.55946 
5409 0.8207 —1.5828 0.051243 —0. 022321 0.1723 33.83 26.31 0.56003 
-6761 0.8087 —0.9865 0.055472 —0.021025 0.08085 31.96 14.07 0.56070 
- 8452 0.8010 —0.6049 0.052354 —0.034415 0.03709 30.79 9.014 0.56114 
014 0.7969 —0.3976 0.051167 —0.032191 0.01930 30.17 7.179 0.56137 
.183 0.7943 —0.2726 0.056409 —0.031203 0.01097 29.79 6.383 0.56151 
.352 0.7927 —0.1914 0.063798 —0.0s7135 0.026700 29.55 5.991 0.56160 
-690 0.7908 —0.09601 0.061576 —0.0s2963 0.022887 29.27 5.655 0.56171 
.028 0.7897 =—0. 04418 0.077653 —0.011439 0.021435 29.12 5.530 0.56177 
0.7874 and a2=0.07361. The b's and ec’s are independent of A. 
Table 21 
Coefficients Kst, Ki, and Kym, (a/h = 1000) 
A=0 =r A=0 
n Ka X 108 Ki X 106 Kat X 108 Ki X 108 Kem X 108 
(Von Mises 
Formula) 
20 63.38 91.66 57.50 85.24 88.66 
19 50.89 68.80 47.05 65.19 
18 47 .83 67.86 44.51 64.29 65 .93 
17 39.33 52.87 37.21 50.05 
16 36.61 52.15 34.64 50.02 50.70 
14 25.71 34.05 24.82 33.14 
13 23.88 33.59 23.06 32.70 
12 22.71 33.95 21.92 33.05 
12 vee 22.37 17.38 21.98 21.77 
11 16.20 21.71 15.85 21.33 
10 15.07 21.76 14.73 21.38 
10 12.34 15.48 12.15 15.29 
9 ale aly 15.08 11.00 14.89 
9 9.59 11.65 9.49 11.54 11.33 
8 8.576 11.130 8.480 11.03 
8 7.493 8.983 7 A27 8.920 
te 6.547 8.532 6.489 8.472 
6 4.466 5.543 4,442 5.517 5.356 
5 3.963 5.569 3.941 5.5438 
5 3.102 3.843 3.089 3.830 
5 2.611 2.882 2.605 2.876 
4 1.967 2.445 1.962 2.439 1.973 
4 1.597 1.715 1.595 1.713 
3 1.094 1.372 1.092 1.370 
3 0.9090 0.9958 0.9080 0.9950 0.9728 
3 0.8372 0.8600 0.8366 0.8595 
2 0.4479 0.6079 0.4475 0.6074 
2 0.3566 0.3795 0.3564 0.3793 0.3722 
2 0.3172 0.3173 0.3172 0.3172 
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Coefficients Kst, Ki, and Kym, (a/h = 100) 


A=0 = A=0 
n Kat X 108 n Ki X 108 n Ks X 105 n K; X 106 Kom X 108 
20 9273 
19 12120 19 8454 
15 7577 .6 
18 11020 18 7684 
14 12000 14 6707 
17 10000 17 6969 
16 9055 16 6310 5344 
14 7432 14 5175 4544 
13 6484 13 4161 
12 5731 12 3651 
13 4813 13 3693 3530 
11 3484 11 2473 
11 2500 11 2134 2185 
10 1905 10 2268 10 1544 10 1935 
8 1402 10 1663 8 1132 10 1498 1548 
9 1481 9 1334 
8 1041 8 904.4 
7 882.0 if 765.1 
8 1039 8 964.0 
th 712.0 i 644.1 
6 605.5 7 753.8 6 547.2 7 712.9 
7 648.5 ui 621.3 634.7 
6 506.0 6 469.5 
6 451.7 6 557.4 6 533.6 
5 425.5 5 394.0 
5 361.5 5 340.9 
6 435.6 6 423.5 
5 297.0 5 360.0 5 286.1 5 349.8 
4 230.4 5 285.2 t 221.6 5 279.9 281.1 
4 186.4 4 226.7 4 181.9 4 222.3 
4 167.1 4 183.6 4 164.4 4 181.1 180.6 
3 140.5 3 136.8 
4 163.9 4 162.3 
3 113.4 3 151.0 3 111.4 3 148.8 
3 98.93 3 116.6 3 97 .62 3 115.3 
3 91.13 3 99.41 3 90.23 3 98.59 97.46 
3 83.28 3 84.57 3 84.13 
2 71.08 2 69.81 
2 56.08 2 55.32 
2 42.96 2 54.62 2 42.58 2 54.23 
2 37 .22 2 40.76 2 37.00 2 40.56 
2 32.89 2 33.02 2 32.79 2 32.91 32.59 
2 29.28 2 29.22 4 


